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Abstract 

The data-random graphs called proximity catch digraphs (PCDs) have been introduced recently and have 
applications in pattern recognition and spatial pattern analysis. A PCD is a random directed graph (i.e., 
digraph) which is constructed from data using the relative positions of the points from various classes. Different 
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the underlying graphs based on a family of PCDs which is determined by a family of parameterized proximity 
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the underlying graphs. We demonstrate that, properly scaled, relative edge density of the underlying graphs 
is a {/-statistic, and hence obtain the asymptotic normality of the relative edge density for data from any 
distribution that satisfies mild regulatory conditions. By detailed probabilistic and geometric calculations, we 
compute the explicit form of the asymptotic normal distribution for uniform data on a bounded region. We 
also compare the relative edge densities of the two types of the underlying graphs and the relative arc density 
of the PCDs. The approach presented here is also valid for data in higher dimensions. 
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1 Introduction 

Classification and clustering have received considerable attention in the probabilistic and statistical literature. 
In this article, the probabilistic properties of a graph invariant of a family of random graphs is investigated. 
Data-random digraphs are directed graphs in which each vertex corresponds to a data point, and directed edges 
(i.e., arcs) are defined in terms of some bivariate relation on the data. For example, nea r est ne ighbor digraphs 
are defined by placing an arc between each vertex and its nearest neighbor. iPriebe et al.l (|200ll ) introduced the 
class cover catch digraphs (CCCDs) in R and gave the exact and the asymptotic distrib u tion of the domination 



numb e r of the CCCDs for uniform data on bou nded intervals. iDeVinnev et al.l (120021) , iMarchette and Priebe 



(|2003l) . IPriebe et al J (|2003ah . IPriebe et al.l (|2003rjt) . and lDeVinnev and Priebel (|2006t applied the concept in higher 
dimensions and demonstrated relatively good performance of CCCDs in classification. Their methods involve data 
reduction (i.e., condensing) by using approximate minimum dominating sets as prototype sets (since finding the 
exact minimum dominating set is an NP-hard problem in general — e.g., for CCCD in multiple dimensions — (see 
DeVinnev and Prie bel (|200fih ). Furthermore the exact and the asymptotic distribution of the domination number 
of the CCCDs are not analytically tractable in multiple dimensions. For the dom ination number of CCCDs for 
one-dimensional data, a SLLN result is proved in IDeVinnev and Wiermanl (120031 ). and this result is extended by 
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Wierman and Xiang <l2008h : furthermore, a generalized SLLN result is provided by IWierman and Xiana (|2008l ). 



and a CLT is also proved bv lXiang and Wierman ( 20091 ). The asymptotic distribut ion of th e dom ination number 
of CCCDs for non-uniform data in R is also calculated in a rather general setting ( Cevhanl ( 20081 )). 

Cevhan (|2005l) generalized CCCDs to what is called proximity catch digraphs (PCDs). Let (fl,A4) be a 
measurable space and X n — {X\, X2, ■ ■ ■ , X n } and y m = {Yi, Y2, . . . , Y m } be two sets of Sl-valued random 
variables from classes X and y, respectively, with joint probability distribution Fx,y- A PCD is comprised by 
a set of vertices and a set of arcs. For example, in the two class case, with classes X and y, the X points are 
the vertices and there is an arc from x\ € X n to X2 € X n , based on a binary relation which measures the relative 
allocation of x\ and X2 w ith respect to y points. The PCDs are closely related to the class cover problem of 
Cannon and Cowen (120001) . The class cover problem for a target class, say X, refers to finding a collection of 
neighborhoods, Ni around X, such that (i) X n C (\J i Ni) and (ii) y m r\(\J i Ni) = 0. A collection of neighborhoods 
satisfying both conditions is called a class cover. A cov er satisfying ( i ) is a prope r cover of class X wh i le a c over 
satisfying (ii) is a pure cover relative to class y. See IPriebe et all (|200ll ) and Cannon and Cowe pi (120001) for 
more detail on the class cover problem. The firs t PCD family i s intr oduced bv lCevhan and Priebd ( 2003 ); the 
parametrized version of this PCD is developed bv lCevhan et al.1 (2007 ) where the relative arc density of the PCD 
is calculated and used for spatial pattern analysis. Cevhan and Priebe ( 20051 ) introduced another digraph family 
called proportional ed ge PCDs and calcula t ed th e asympto t ic dis tribution of its domination number and used it 
for the same purpose (jCevhan and Priebd (|20071 ); ICevhan ( 201C)). The re lative arc density of this PCD family 



is also computed and used in spatial pattern analysis (jCevhan et al. (2006)) 



The underlying graphs based on digraphs are obtained by replacing arcs in the digraph by edges based on 
bivariate relations. If symmetric arcs are replaced by edges, then we obtain the AND-underlying graph; and if 
all arcs are replaced by edges without allowing multi-edges, then we obtain the OR-undcrlying graph. Properly 
scaled, we demonstrate that the relative edge density of the underlying graphs of proportional-edge PCDs is a 
[/-statistic, which has asymptotic normality by the general central limit theory of [/-statistics. Furthermore, we 
derive the explicit form of the asymptotic normal distribution of the relative edge density of the PCDs based on 
uniform data in a bounded region. 

For the digraphs introduced by IPriebe et all (|200ll ) (i.e., CCCDs), whose relative arc density is also of the 
U -statistic form, the asymptotic mean and variance of the relative den sity is not ana l ytical ly tra ctable, due to 
geome tric difficulties encountered. However, for the PCDs introduced in ICevhan et all (|2006f ). and[c evhan et alJ 
( 2007 ). the relative arc density has tractable asymptotic mean and variance. We define the relative densities of 
graphs and digraphs and derive their asymptotic distribution in general in Section [2] define the underlying graphs 
of proportional-edge PCDs and their relative edge densities in Section [3l provide the asymptotic distribution of 
the relative edge density for uniform data in Section [4] We treat the multiple triangle case in Section [5l provide 
the discussion and conclusions in Section [6l and the tedious calculations and long proofs are deferred to the 
Appendix. 



2 Relative Density of Graphs and Digraphs 

The main difference between a graph and a digraph is that edges are directed in digraphs, hence are called arcs. 
So the arcs are denoted as ordered pairs while edges arc not. 

2.1 Relative Edge Density of Graphs 

Let G n = (V, £ ) be a graph with vertex set V = {v\, V2, ■ • ■ , v n } and edge set £. The relative edge density of the 
graph G n which is of order |V| = n, denoted p(G n ), is defined as 

P{G n ) 



n(n — 1) 



where | • | denotes the set cardinality function ( Janson et al.1 ( 2000| )). Thus p{G n ) represents the ratio of the 



number of edges in the graph G n to the number of edges in the complete graph of order n, which is n(n — l)/2. 
If G n is a random graph in which edges result from a random process, the edge probability between vertices Vi, Vj 
is defined as p e {i, j) '■= P(viVj € £) for all i 7^ j, i,j = l,2,...,n. 

Theorem 2.1. (Main Result 1) Let G n = (V,£) be a graph of order n with V = {vi, V2, ■ ■ • , v n } and let 

ht 3 := Hv iVj S B ). 



(a) If the set £ of edges result from a random process, then p{G n ) is a one-sample U -statistic of degree 2. 
Moreover, if p e (i, j) = p e for all i ^ j, i, j = 1, 2, . . . , n (i.e., the edge probability is constant for each pair 
of vertices Vi,Vj), then p{G n ) is an unbiased estimator ofp e . 

(b) If the set £ of edges result from a random process, such that hfj are identically distributed with p e (i,j) = p e 
for alii jtz j, i,j = 1,2, ... ,n, and hfj and h kl are independent for distinct i,j,k,l, and v e := Cov(^, hf k ) > 

for all i =/= j =/= k, i,j, k = 1, 2, . . . ,n, then y/n [p(G n ) — p e ] jV(0, 4 v e ) as n — > oo, where — stands 
for convergence in law or distribution and J\f(p,a 2 ) stands for the normal distribution with mean p and 
variance a 2 . 

Proof: (a) Assume the edges £ result from a random process and let G n be the corresponding graph. Let 
hfj = l(viVj £ £). Since the edge ViVj £ £ can equivalently be expressed as VjVi £ £ for all we have hfj = h e j i 
and so hfj is symmetric in Additionally, \£\ — J2i<j hfj- So 



p{G n ) — tttt Kj- 

Thus. p{G n ) is a one-sample {/-statistic of degree 2 with symmetric kernel hfj. Assume, moreover, P(viVj £ £) = 
p e for all i ^ j, i,j = 1,2, ... , n. Then for i ^ j, we have E[hf 3 ] = E[/if 2 ] = E[I(viu 2 £ £)] = P{viv 2 £ £) = p e 
Hence p e is an estimable parameter of degree 2. Furthermore, 

i<j i<j 



E[p(G„)] = -j^— E[|£|] = £ m>] =Ef- -T^S ^T^Pe = Pe- 
n{n— 1) f . n(n — 1) V 2 / 



Then, p(G n ) is actually an unbiased estimator of p e . 

(b) Assume the conditions for G n = (V,£) stated in the hypothesis. In part (a) we have shown that p e is an 
estimable parameter of degree 2, and p{G n ) is a one-sample [/-statistic of degree 2 with symmetric kernel hfj . Fur- 
thermore, (hfj) 2 = hfj, since {I{v tVj £ £)) 2 = I(viVj £ £). So E[(^) 2 ] = E[ft?-] = p e < oo and Cov^-, hf k ) = 
EKMk] ~pI< o o, since E^-feffc] = P((hf 3 ,hf k ) = (1, 1)). By the hypothesis, v e = CawQifj, hf k ) > 0. Then by 

Theorem 3.3.13 in lRandles and Wolfel (|l979l ). we have y/n [p{G n ) - p e ] J\f(Q,4v e ) asn^oo. ■ 

In part (b) of Theorem O we have Cov(hfj , hf k ) = V[hfjhf k \ - E[^]E[^ fe ] = E[hfjhf k ] ~ p 2 e where 
h ij h ik = S £)I(viVk £ £) = WiUfe} C £), so E[/i^7i^ fe ] = P^u,-, UjUfc} C £). Hence ^ e > iff 

P{{ v i v ji v i v k} C £) > p 2 . Notice that E[|/if-| 3 ] = E\hfj] = p K < oo and assuming v P > 0, then the sharpest rate 



of convergence in the asymptotic normality of p(G n ) is ( Callaert and Janssen (Il978h ) as follows 



sup 



p (Vn-(p(G n )- Pe ) <f] _ m 
V 2^ 



<C .p e •(4^ e )- 3 / 2 -n- 1 /2 = C e 



where C e is a constant and $(•) is the standard normal distribution function. Furthermore, we have 
Var[ft?.] = E[(^ 7 ) 2 ] - (E[/if-]) 2 - E[ft?-] -p 2 = p e - p\ = p e (l - p e ). 



The graph G n in Theorem 12. II is not a deterministic graph, but a random one. In general a random graph is 
obtained by starting with a set of n vertices and adding edges between them at random. Most commonly studied 
is the Erdos-Renyi model, denoted G(n,p), in which every possible edge occurs independently with probability 
p ( Erdos and Renvi (|l959f )). Notice that the random graph G(n,p) satisfies part (a) of Theorem 12.11 so the 



relative edge density of G(n,p) is a [/-statistic; however, the asymptotic distribution of its relative edge density is 

degencrati 
the edges. 



degenerate (with p(G(n,p)) — > pasu-) oo) since the covariance term is zero due to the independence between 



2.2 Relative Arc Density of Digraphs 

Let D n = (V, .4) be a digraph with vertex set V = {v±, v 2 , ■ ■ ■ , v n } and arc set A. The relative arc density of the 
digraph D n which is of order |V| = n, denoted p(D n ), is defined as 



Thus p(D n ) represents the ratio of the number of arcs in the digraph D n to the number of arcs in the complete 
digraph of order n, which is n(n — 1). If D n = (V,A) is a random digraph in which arcs result from a random 
process, the arc probability between vertices Vi,Vj is defined as p a (hj) '■= P(( v ii v j) € A) for all i 7^ j, i,j = 
1,2,. ..,n. 

Theorem 2.2. (Main Result 2) Let D n = {V,A) be a digraph of order n with V = {v\, v-x, ■ ■ ■ , v n } and let 
gij := I((vi,Vj) G A). 

(a) If the set A of arcs result from a random process, then p{D n ) is a one-sample U -statistic of degree 2. 
Moreover, if p a (i,j) = Pa for all i 7^ j , i,j = 1,2, ... ,n, (i.e., the arc probability is constant for each pair 
of vertices Vi,Vj), then p(D n ) is an unbiased estimator of p a - 

(b) If the set A of arcs result from a random process such that g^ are identically distributed with p a (i, j) = Pa for 
all i ^ j, i,j = 1,2, ... , n, g^ and g k i are independent for i ^ k and j 7^ I (i.e., gij and gki are independent 
for distinct i,j,k,l and gij and gji are independent), and Cov{gij, gu) > for all i 7^ j and k 7^ I and 

exactly one of i, j is equal to exactly one of k, I for i, j, k, I = 1, 2, . . . , n, then ^fn [p(D n ) — p a ] — > A/"(0, v a ) 
as n — > 00, where v a = lim,,-^ n Var [/>(£>„)] . 

Proof: (a) Assume that the arcs A result from a random process and let D n be the corresponding digraph. Let 
gij = I((vi,Vj) G A). The arcs (vi, Vj) £ A and (vj, Vi) G A are distinct for i 7^ j, so g^ is not symmetric in 
But we can define a symmetric kernel as hfj = (g^ + gji)/2. Then we have, — J2i<j ^tj- So 

P( D n) = ?4tX! 

{ 2 > i<i 

Thus, p(D n ) is a one-sample [/-statistic of degree 2 with symmetric kernel hfj. Assume, moreover, P((vi,Vj) G 
A) = p a for all i ^ j, i,j = 1,2, ... , n. Then for i 7^ j, we have 

E[/i£] = E[( 5 , y + gji )/2) = (E[ 9ij ] + E[fcd)/2 = Efoy] = E[ 5l2 ] = P((v 1} v 2 ) G A) = Pa . 

Hence p a is an estimable parameter of degree 2. Furthermore, 

n[n — 1) n(n — 1) ^— ' J n(n — 1) z —' n(n — 1) ' 

i<3 i<j i<j 

Then, p(D n ) is actually an unbiased estimator of p a . 

(b) Assume the conditions for D n = (V,A) stated in the hypothesis. In part (a) we have shown that p a is an 
estimable parameter of degree 2, and p{D n ) is a one-sample [/-statistic of degree 2 with symmetric kernel hfj. 
Furthermore, (h^) 2 = { g 2 J +2g lJ gj l +g 2 t )/A = (gij+2g ij g ji +g ji )/4, since (I((vi,Vj) G A)) 2 = I((vi,Vj) G A). So 
E[(^) 2 ] = (Efotf] + 2 K //, ,//,, + = (2 Pa + 2 E[ff <iflii ])/4. Since E[s y5ii ] = E[I((^-) G A)I(( Vj , v t ) G 

.A)] = E[I({(u l ,w J ), (wj,^)} C A)] = P({(vi,Vj), (vj,Vi)} C A), we have E[(^) 2 ] < 00. By the hypothesis, v a = 
Cov(hfj,hf k ) = Cov((g ij +g j i)/2,(g i k+gki)/2) = {Cow[g l3 , g tk ]+Cow[g lJ , g M ]+Cow[g Jt , g lk ] + Co\[g Jt , g kl ]) / A > 
0. Additionally, Cov( gij ,g k i) = ^[g^gu] - p 2 a < 00, since E[g>yg w ] = P{{g i3 ,g k i) = (1,1))- Hence ^ a < 00 as 
well. Then by Theorem 3.3.13 in iRandles and Wolfd (Il979l) . we have ^/n [p{D n ) —p a ] — > A/"(0, v a ) as n — > 00. ■ 

In part (b) of Theorem 12. 2[ we have 



Cov(^-, /&) - E[ft?.fc? fc ] - E[^-]EK fc ] - E[/i?./ifj.] - pi 



where 



K 3 K k = (9ij + 9]i){9ik + 9ki)/4: = (gijg ik + gijg ki + g^ga, + gji9ki)/^ = (!{{vi,Vj) G A)I((vi,v k ) G A)+ 
!{{Vi,Vj) G A)I{{v k ,Vi) G A) + I((u J -,u i ) G A)I((ui,v fe ) G A)+I{{vj,Vi) e A)I((v k , Vi ) G -4))/4 = 

(UXO^U;), (Ui,Vfe)} C A)+I({(w l ,Wj):K,Wi)} C A)+I({(Vj,V l ),(v l ,V k )} C A)+I({(Wj,U i ),(Wfe,Wi)} C A))/A. 

So, ^ a > iff 

(^({(ui.^J.fwi.Wfc)} C A) + P({(v i ,v j ),{v k ,Vi)} C ^)+P({(w J -,« i ),(«i,«*)} C ^) + P({fe,^),K,^)} C A))/4>p 2 a . 
Notice that 

E[|^| 3 ] - E[( ftj + g ri ) 3 /8] - E[4 +34^ • 3 + 5 | 4 ]/8 = +3^-^ + 3g ijgji + g Tl }/8 = 

(2E[ff tf ] + 6E[3y 5j i])/8 = (pa + 3p SQ )/4 < 00, 



where p sa = P(gijfjji = 1) = P({(uj, Vj), (vj,Vi)} C A) is the symmetric arc pro bability in D n . Assum i ng v a > 0, 
then the sharpest rate of convergence in the asymptotic normality of p(D n ) is ( Callaert and Janssen ( 19781 )) as 
follows: 

Pa 



n ) Pa ) 



sup 

tGK 



where C a is a constant. Furthermore, we have 



-*(*) 



< C ■ Pa ■ Kr 3/2 • n- 1 / 2 = C a 



Var[/£] = E[(h? 



E[.9, 



%jj j 

2 9ij9ji 



"5ii]/4- 



E[(Py+5ii)/4] 



E k 2 , + 2gij9ji 



pI 



(Eb«] + 2E[ ff „ Sji ]+E[ Sji ])/4- 



"4l/ 4 - Pa 



■p«)/4- 

{Pa+Psa)/2 



pI 



The digraph Z? n in Theorem 12.21 is not a deterministic digraph, but a random one. In general a random 
digraph, just like a random graph, can be obtained by starting with a set of n vertices and adding arcs between 
them at random. We can consider the counterpart of the Erdos-R enyi model for digraphs , denoted D(n,p), 
in which every possible arc occurs independently with probability p ( Erdos and Renvi ( 19591 ) ). Notice that the 
random digraph D(n,p) satisfies part (a) of Theorem 12. 2[ so the relative arc density of D(n,p) is a {/-statistic, 

however, the asymptotic distribution of its relative arc density is degenerate (with p(D(n,p)) — > p as n 
since the covariance term is zero due to the independence between the arcs. 



oo 



3 Relative Edge Density of Underlying Graphs Based on PCDs 

3.1 Proximity Catch Digraphs and the Corresponding Underlying Graphs 

Let be a measurable space and d(-,-) : O x fi — > [0, oo) be any distance function. Consider N : fl — >• 

"P(f2), where V(-) represents the power set functional. Then given y m C fl, the proximity map N(-) associates 
with each point x £ f2 a proximity region N(x) C ft. The region N(x) is defined in terms of the distance 
between x and y m . Define the data-random PCD, D n , with vertex set V = {X\,X2, ■ ■ ■ ,X n } and arc set A 
by (Xi,Xj) e A Xj 6 N(Xi) where point Xi "catches" point Xj. The random digraph D n depends on 

the (joint) distribution of the Xi and on the map iV(-). The adjective proximity — for the catch digraph D n 
and for the map N(-) — comes from thinking of the r egion N(x) as representing those points in ft "close" to x 
( Toussaintl ( 1980h and ljaromczvk and Toussaint ( 19921 )). The Ti-region Ti(-, N) : Q — s- V(fl) associates the region 



Ti(x,N) := {z e O : x £ N(z)} with each point x £ Q. A Ti-rcgion is sort of a "dual" of the corresponding 
proximity region and is closely associated with domination number being equal to one. If X\, X2, ■ ■ ■ ,X n are 
^-valued random variables, then the N(Xi) (and T±(Xi, N)), i = 1,2, ...,n are random sets. If the Xi are 
independent and identically distributed, then so are the random sets N(Xi) (and r±(Xi, N)). 

If Xi, X2, ■ • ■ , X n ~ F, then, by Theorem l2.21 the rel ative arc density of the asso ciated data-random proximity 



catch digraph, D n , denoted p(D n ), is [/-statistic. See ICevhan et al.l (|2007l 120061 ) for its derivation and other 
details. 

The underlying graph of a digraph is the graph obtained by replacing each arc (it, v) £ A or each symmetric 
arc, {(«, v), (v, u)} C A by the edge uv. The former underlying graph will be referred as the OR-underlying graph, 
while the latter as the AND-underlying graph. That is, the AND-undcrlying graph for digraph D n = (V,A) is 
the graph G an d(D n ) = (V,£ a nd) where £ an d is the set of edges such that uv £ £ anc j iff (u,v) £ A and (v,u) £ A. 
The OR-underlying graph for D n = (V, A) is the graph G OI (D n ) = (V, £ OT ) where £ OI is the set of edges such that 
uv £ for iff (it, v) £ A or (v, u) £ A. 

Consider the data-random PCD, D n , with vertex set V = {X\, X2, ■ ■ . , X n } and arc set A defined by (Xi, Xj) £ 
A <i=J> Xj £ N(Xi). The AND-underlying graph, G an d, of D n with the vertex set V and the edge set £ an( j 
is defined by XiXj £ £ anc i iff (Xi,Xj) £ A and (Xj,Xi) £ A. Likewise, the OR-underlying graph, G or , of D n 
with the vertex set V and the edge set £ or is defined by XiXj £ £ OI (Xi,Xj) £ ^4 or (Xj,Xj) £ ^4. Then 

JQX; £ f and iff £ N(Xi) and X, £ N{Xj) iff £ X(A 4 ) and Xj £ Li(AT t , X) iff Xj £ iV(Xi) n Li(X,, N). 
Similarly, X,X 3 £ £ or iff Xj £ iV(Xi) U T^X,, N). 



3.2 Relative Arc Density of the Proximity Catch Digraphs 



The relative arc density of the PCD, D n , is denoted as p(D n ). Let g^ = I((Xi,Xj) € A) — I(Xj <G 7V(Xj)). 



Then for Xj ~ F, i = 1, 2, . . . , n, p(D n ) can be written as 

2 



p(D n ) 



n (n — 1) 



i<j 



where 2 hij = (gij + gji) = I(Xj E N(Xi)) + I(X,; <G N(Xj)) is the number of arcs between X 2 ; and Xj in D n . 
Note that hij is a symmetric kernel with finite variance since < hij < 1. Moreover, p(D n ) is a random variable 
that depends on n, F, and N(-) (i.e., 3^)- But E \p(D n )] only depends on F and N(-). That is, 



< E [p(D n )] = — £ ]T Elhij] = E [h 12 ] 



(1) 



where 2E[ft 12 ] = E[I((X!,X 2 ) G .A) +I((X a ,Jri) G A)} = E[I((X U X 2 ) G A)} +E[I{{X 2 ,X 1 ) G 4)] = 
P((Xi,X 2 ) G .A) + P((X 2 ,X 1 ) G A) = P(X 2 G Ar(Xi))+P(Xi G AT(X 2 )) = 2 Po (Af)- Hence E [hi 2 ] = p a (N), 
which is the arc probability for the PCD, D n . Notice also that p a (N) = P(Xj G N(Xi)) for i ^ j. Furthermore, 



< Var [p(D n )} 



n 2 (n — 1)' 



■Var 



EE* 



i<3 



(2) 



Expanding this expression, we have 



2 4(n — 2) 

Var [p(D n )\ = — -Var [h 12 ] + — ± Cov [h 12 , h 13 ] . 

n(n — l) n (n — 1) 



As in Section l2~2l we have 



Var [h 12 ] = ( Pa (N) ~ [ Pa (N)} 2 ) = Pa (N) (1 - Po (JV)) ■ 

Moreover, the covariance is as follows 

Cov [h l2 , h 13 ] = E [h 12 .h 13 ] ~ E [ha] E [h 13 ] , 

where E [hi 2 ] — E [hi 3 ] — p a (N) and, 

4E [hi 2 .hi 3 ] = E[(gi 2 + g 2 \)(gi 3 + g 3 i)} = E[gi 2 gi 3 + .912531 + 521.913 + .921.931] 

= E[I(X 2 G N{Xi)l{X 3 G N(Xi)) + I(X 2 G N(Xi)I(Xi G JV(X 3 )) + 

I(Xi G iV(X 2 )I(X 3 G JV(JCi))] + G iV(X 2 )I(Xx G AT(X 3 ))] 
= E[I({X 2 ,X 3 } C iV(Xi)) + I(X 2 G iV(X!)I(X 3 G ri(jr 3 ,JV)) + 

I(X 2 G ri(Xi)I(X 3 G N(X 1 )))+I(X 2 G r^Xr.^ICXa G ^(X^iV))] 
= P({X 2 ,X 3 } c JV(Xi)) +2P(X 2 G 7V(Xi),I(X 3 G ri(X 3 , TV)) + P({X 2 ,X 3 } C ri(X!,iV)). 

The digraph D n is a random digraph where the arc probability is P((Xj,X 7 ) G A.) = p a (N) for i 7^ j and 
is an estimable parameter of degree 2. Using Equation (JXJ) , we have that p(D n ) is an unbiased estimator of 
p a (N). Notice that for PCDs, the set of vertices V = X n is a random sample from a distribution F (i.e., the 
vertices directly result from a random process), and the arcs arc defined based on the random sets (i.e., proximity 
regions) N(Xi) as described before. Hence the set of arcs A (indirectly) result from a random process such 
that g^ are identically distributed and g^ and gki are independent for i ^ k and j ^ I. Furthermore, we have 
v a {N) := Cov [hij, hik] < 00 as before. Then we have the following corollary to the Main Result 2. 

Corollary 3.1. The relative arc density, p(D n ), of the PCD, D n , is a one-sample U -statistic of degree 2 and is 
an unbiased estimator of p a (N). If, additionally, v a (N) = Cov [hij, h^] > for all i 7^ j =/= k, i,j, k = 1, 2, . . . , n, 

then sfn [p(D n ) - p(N)} -A Af(0, 4 v a (N)) as n -> 00. 



In the above corollary, v a (N) > iff P({X 2 ,X 3 } C iV(Xi)) + 2P(X 2 G iV(Xi),I(X 3 G Ti(X 3 ,N)) + 

p({x 2 ,x 3 } c ri(x lf iv)) > 4 P 2 a (N). 



3.2.1 The Joint Distribution of 

The pair (/112, his) is a bivariate discrete random variable with nine possible values: 

(hi2,h 13 )€ {(0,0), (0,1), (0,2), (1,0), (1,1), (1,2), (2,0), (2,1), (2, 2)}. 

Then finding the joint distribution of (ft.12,^13) is equivalent to finding the joint probability mass function of 
(h 12 ,hi 3 ). 

First, note that 

(hi2, hi 3 ) = (0, 0) iff gi2 = 321 = .913 = .931 = iff 

l(x 2 g N(Xi)) = I(Xi e N(x 2 )) = I(X 3 e JV(*i)) = l(Xi e N(x 3 )) = iff 
i(x 2 g n(Xi)) = i(x 2 g r^XuN)) = i(x 3 g jv(Xi)) = i(x 3 g r^, AO) = o iff 
i(x 2 g r(y 3 )\A r (^i)) = i(x 2 e T(y 3 )\ri(X!,iV)) = I(X 3 G T{y 3 )\N(X x )) - I(X 3 G T(y 3 )\r 1 (x 1 ,7V)) = 1 iff 
i({* 2 ,* 3 } c T(y 3 )\N(x 1 )) = i({x 2 ,x 3 } c T(y 3 ) \ t^x^n)) = 1 iff 
i({^ 2 ,^ 3 } c T(y 3 ) \ (Nix,) u iyx^iv))) = 1. 

Hence P((h 12 ,h 13 ) = (0,0)) = P({X 2 ,X 3 } C T(y 3 ) \ (JV(Xi) U r^Xi, AO)). 

Furthermore, by symmetry, P((h 12 ,h 13 ) = (0, 1)) = P((hi2, his) = (1, 0)), P((/na, h 13 ) = (0, 2)) = P((h 12 ,h 13 ) = 
(2, 0)), and P((hi 2 , hi 3 ) = (1, 2)) = P((hi 2 , h\ 3 ) = (2, 1)). So it suffices to calculate one of each pair of the prob- 
abilities in the above cases. 

Finally, 

(hi 2 , hi 3 ) = (2, 2) iff gi2 = .921 = .913 = .931 = 1 iff 
I(X 2 G N(X X )) = I(X a G AT(X 2 )) = I(X 3 G iV(Xi)) = I(Xi G X(X 3 )) = 1 iff 

i(x 2 g iv(Xi)) = i(x 2 g r x (Xi,iv)) = i(x 3 g iV(Xi)) = i(x 3 g t^x^n)) = 1 iff 
i({x 2 ,x 3 } c iv(xo) = i({x 2 ,x 3 } c r x (xi,jv)) = 1 iff 
i({x 2 , x 3 } c (jv(Xi) n ri(Xi, AO)) = 1. 

Hence P((h 12 ,h 13 ) = (2,2)) = P{{X 2l X 3 } C [N(X l )^T 1 {X 1 ,N))). Finally, P((h 12 ,h 13 ) = (1,1)) can be found 
by subtracting the sum of the probabilities in the other cases from 1. 

3.3 Relative Edge Density of the AND-Underlying Graphs Based on PCDs 

The relative edge density of the AND-underlying graph, G a nd(-Dn), based on the PCD, D n , is denoted as Pa, n d(D n ). 
For Xv ~ F, i = 1,2,... , n, one can write down the relative edge density as 

Panel (A*) = -, 2 - V V 

where 

h™ d = I((Xi,Xj) G f and ) = I((X, X,-) G A, {X jt Xi) G A) 

= I(Xj G JV(Xi), X G N(Xj)) = I(Xj G A^(X), X,- G Ti (X», AT)) 

= i(Xj eiV(x)nri(x,A^)) 

is the number of edges between Xi and Xj in G an d(-Dn) or number of symmetric arcs between Xj and Xj in 
Note that /ifj ld is a symmetric kernel with finite variance since < ft.f" d < 1. Moreover, p a nd(-Dn) is a random 
variable that depends on n, F, and N(-) (i.e., y). But E [p an d(-Dn)] only depends on F and AT(-). That is, 

< E [ Pand (D n )] = - 2 _ 2 2 E[^ d ] = E [/^ d ] (3) 

^ ' i<3 

where E [hf^\ = E[I((Xi,X 2 ) G A , (X 2 ,X 1 ) G 4)] = P((X lf X 2 ) G .4 , (X 2 ,Xi) G .4.) = P(X 2 G X(X X ) n 
Fi(Xi, AT)) = p ail< i(N). Notice that p an d(N) is the edge probability for the underlying graph G an <i(D n ), but it is 



symmetric arc probability p sa {D n ) for the PCD, D n . Notice also that Pand(-W) = P(Xj G N(Xi) n Ti (Xi,N)) 
for i =/= j- Furthermore, 



< Var [p and (D n )} = 4 ^ 2 Var 
n z (n — l) z 



E E ^ 



i<3 



(4) 



Expanding this expression, we have 



Var [ Pand (A0] = ^TT) Var [^» d ] + i^-^ Cov [h™ d , h% d ] . 

Here, as in Section 12.11 we have 

Var [hlf] = Pand (N) - [p ancl (iV)] 2 = p and (iV) (1 - Pand (JV)) . 
Moreover, the covariance is as follows 

Cov [htf, Kf] = E [Kf.hlf] - E [h™ d ] E [h^f] • 
Since E [h^f] = E [hfi*] = p aBd (N) and, 

e [h^ d .hif] = E[i(x 2 G JV(x x ) n^XuN)).!^ e jv(Xi) nri(Xi.JV))] 
= E[i(x 2 € iVpronriOYi.-ZV), x 3 e iV(Xx)nri(Xx,iV))] 
= p(x 2 g ^(Xijnri^.JV), x 3 g AT(Xi)nr 1 (x 1 ,iV)) 

= P({X 2 ,X s }cN(X l )nT 1 (X l ,N)) 



it follows that 



Cov [htf, h$f] = p({x 2 ,x 3 } c NiXt) n r^Xx, jv)) - bwdWf 



The underlying graph G an d(-Dn) is a random graph where the edge probability is P(XiXj G £ an d) = Pand(N) 
for i 7^ j and is an estimable parameter of degree 2. Using Equation we have that G an d(D n ) is an unbiased 
estimator of Pa, n d(N). Notice that for the AND-underlying graphs based on the PCDs, the set of vertices V = X n is 
a random sample from a distribution F (i.e., the vertices directly result from a random process), and the edges are 
defined based on the random sets N(Xi) H Ti (Xi,N) as described before. Hence the set of edges £ an d (indirectly) 
result from a random process such that /ifj ld are identically distributed and hf^ d and hf.f d are independent for 
distinct k, I. Furthermore, we have Cov [hfj d , /i^" d ] < oo as before. Then we have the following corollary to 
the Main Result 1. 

Corollary 3.2. The relative edge density, p an d{Dn)j of the AND-underlying graph, G an d(D n ), is a one-sample 
U -statistic of degree 2 and is an unbiased estimator of p an( i(N) . If, additionally, v a nd{N) '■= Cov [h° , J ld , h° , j^ d ^ > 

for all i^ j ^ k, k = 1,2, ...,n, then sfn [p an d{D n ) - p an d{N)] —> Af(0,Av and (N)) as n ->• oo. 
In the above corollary, iw(iV) > iff P({X 2 ,X 3 } C 2V(X X ) fl T^X^N)) > [ Pa nd(N)} 2 . 

3.3.1 The Joint Distribution of (hf% d , h%£ d ) 

By definition (hf 2 d , ^x3 d ) i s a discrete random variable with four possible values: 

(hlf,hlf) e {(0,0), (1,0), (0,1), (1,1)}. 

Then finding the joint distribution of (h^ , h^ d ^j is equivalent to finding the joint probability mass function of 

( Land ^EindA 
\ n X2 ' 13 )■ 

First, note that 

(ftxTUxf) = (o,o) iff 

\[X 2 € N{X 1 )nT 1 (X 1 ,N)) = I(X 3 G N(Xi) n Tx (X±, N)) = iff 

i(^2 g r(%) \ iV(Xx) n rx(Xx,iV)) = i(x 3 g r(y 3 ) \ JV(*i) n r^x^N)) = 1 iff 

I({X 2 ,X 3 } c T(y 3 ) \ [JV(Xi) n TriX^N)}) = 1. 
Hence P((^ d ,^ d ) - (0,0)) = P({X 2 ,X 3 } C T(y 3 ) \ [JV(X X ) n ^(Xi.iV)]). 



Next, the pair (h^ d , hf^ d ) = (1,1) iff hfi d = hff = 1. So P ((h^ d , hff) = (1, 1)) = E [h^ d .h^ d ] . 
Furthermore, by symmetry P ((h^ d , /ifg d ) = (0, 1)) = P((h^ d , h^ d ) = (1,0)). So it follows that 

P ((^ d , Kf) = (0, 1)) = P ((hlf, hlf) = (1, 0)) 

= \[1-{P {(hlf, hlf) = (0,0)) + P ((htf, h\f) = (1, 1)))] . 

3.4 Relative Edge Density of OR-Underlying Graphs Based on PCDs 

The relative edge density of the OR-underlying graph, G OI (D n ), based on the PCD, D n , is denoted as p m (D n ). 
For Xi ~ F, i = 1, 2, . . . , n, one can write down the relative edge density as 

Por(D n )= 2 h% 
n(n — 1 ' ' J 

where 

= i((x t ,x 3 ) g ^r) = i((X, x,) g .A or (x 3 , g .A) 

= I(Xj G iV(X;) or X; G iV(Xj)) = I(Xj G iV(Xi) or X 3 G Ti (Xi,N)) 
= I(Xj G N(Xi) U i\ (Xj, iV)) 

is the number of edges between Xj and Xj in G or (£>„). Note that /i° r is a symmetric kernel with finite variance 
since < h°* < 1. Moreover, p OI (D n ) is a random variable that depends on n, F, and X(-) (i.e., y). But 
E[/o or (Z? n )] does only depend on F and iV(-). That is, 

< E [p OT (D n )} = £ E E ^°I] = E t^Sl ( 5 ) 

where E[/ig] = E[I(X 2 G iV(Xi) U Ti (Xi,iV))] = P(X 2 G iV(X x ) U T^Xi, N)) = p ot (N). Notice that p or {N) is 
the edge probability for the underlying graph G OT (D n ) and that p or (X) = P(Xj G iV(Xj) U Ti (X it N)) for i 7^ j. 
Similar to the AND-underlying case, we have, 

9 4 (n — 2) 

< Var [p OI (D n )} = — -Var + -) -( Cav[h&, h? a ]. (6) 

n (n — 1) n [n — 1) 

As before, it follows that 

Var [h%\ = p OI (N) - \p OI (N)} 2 = p OI (N) (1 - p OI (N)) . 

and 

Cov [h%, hH] = E [hT 2 .hT 3 ) E {h™} E [hT 3 ] ■ 
Since E [/i° 2 ] = E [hf 3 ] = p OI (N) and, 

^K 2 .hT 3 } = E[i(x 2 GiV(x 1 )ur 1 (Xi,iV)).i(X3G^v(x 1 )uri(x 1) Ar))] 

= E[I(X 2 G X(Xx) U ri(Xi, N) , X 3 G iV(Xx) U ripfi, AT))] 

= p(x 2 g iv(jfi)uri(jfi,jv), x 3 g jv(Xi)uri(Xi,iv)) 
= p({x 2 ,x 3 }civ(x 1 )ur 1 (x 1 ,7V)), 

it follows that 

Cov [fcf a , fcg] = P({X 2 ,X 3 } C N(X t ) U ^(Xx.iV)) - [por(iV)] 2 . 

The underlying graph G ar (D n ) is a random graph where the edge probability is P(XjXj G £ r) = Por(N) for 
i 7^ j and is an estimable parameter of degree 2. Using Equation ([5]), we obtain that G or (D n ) is an unbiased 
estimator of p m (N). Notice that for the OR-undcrlying graphs based on the PCDs, the set of vertices V = X n is a 
random sample from a distribution F (i.e., the vertices directly result from a random process), and the edges are 
defined based on the random sets iV(Xj) U T\{X\,N) as described before. Hence the set of edges £ or (indirectly) 
result from a random process such that h°j are identically distributed and h°j and are independent for distinct 

k, I. Furthermore, we have Cov [ft,° r , < 00 as before. Then we have the following corollary to the Main 
Result 1. 



Corollary 3.3. The relative edge density, p or {D n ), of the OR-underlying graph, G or (D n ), is a one-sample U- 
statistic of degree 2 and is an unbiased estimator ofp or (N). If, additionally, v or (N) :— Cov \h°J, h°£\ > for all 

i 7^ j 7^ k, i, j,k = 1,2, ... ,n, then y/n \p or {D n ) — p or (N)] 7V(0, 4 v or {N)) asn^oo. 
In the above corollary, u OT {N) > iff P({X 2 ,X 3 } C AT(Xx) U T^X^N)) > \p OI (N)] 2 . 

3.4.1 The Joint Distribution of (hfZ,h%£) 

Finding the joint distribution of (h° 2 , h° 3 ) is equivalent to finding the joint probability mass function of (h° 2 , h° 3 ), 
i.e., finding P{{hf 2 , hf 3 ) = for each G {(0, 0), (1, 0), (0, 1), (1, 1)}. 

First, note that 

(hf 2 ,hf 3 ) = (0,0) iff 
I(X 2 G N(X! )L>T 1 (X 1 ,N)) = I{X 3 G N(Xi ) U Tt (X u N)) = iff 

ipG g r(y 3 ) \ ^v(^i) u rx (Xi, n)) = i(x 3 g r(y 3 ) \ u t^x^n)) = 1 iff 
i({x 2 ,x 3 } c T(y 3 ) \ [jv(Xi) u ri(Xi, jv)]) = 1. 

Hence P((^,/if 3 ) = (0,0)) = P({X 2 ,X 3 } C T(y 3 ) \ U ^(Xx,^)]). 

Next, note that (hf 2 , H&) = (1, 1) iff hf 2 = hf 3 = 1. P((h° 12 ,h°l) = (1, f )) = E [^./if 3 ]. 

By symmetry P((h<&,hf 3 ) = (0,1)) = P((^^5s) = (1,0)). Hence 

P((hl 2 ,h? 3 ) = (0,1)) = P(K 2 ,^) = (1,0)) 

= 5 (1 - [iW 2 , ftf 3 ) - (0, 0)) + P((hH,h? 3 ) = (1, 1))]) . 

Remark 3.4. Note that 2^- = hff d + h%, since if g tj = g jt = 0, then 2h tJ = 0, and ft,fj ld = /i° r = 0; if 
9ij = = 1, then 2 ft,^ = 2, and /ifj ld = h°j = 1; and if <?y = and gj, = 1, then 2 /iy = 1, and /if" d = and 
h°j = 1; by symmetry, the same holds when <?jj = 1 and gji = 0. □ 

3.5 Proportional-Edge Proximity Maps and the Associated Regions 

Let ft = M 2 and 3^3 = {yi, Y2, Y3} C M 2 be three non-collinear points. Denote by T(y 3 ) the triangle (including the 
interior) formed by these three points. For r G [l,oo] define N PE {x) to be the proportional- edge proximity map 
with parameter r and T\{x) := Tx (x,N PE ) to be the corresponding Tx-region as follows; see also Figures [1] and 
[2] Let "vertex regions" -R(yx), R(y2), R{y3) partition T(y 3 ) using segments from the center of mass of T(3^ 3 ) to 
the edge midpoints. For x G T(y 3 ) \ y 3 , let v(x) G y 3 be the vertex whose region contains x; x G R(v(x)). If a; 
falls on the boundary of two vertex regions, or at the center of mass, we assign v(x) arbitrarily. Let e(x) be the 
edge of T(y 3 ) opposite v(x). Let t(v(x),x) be the line parallel to e(x) through x. Let d(v(x),£(v(x),x)) be the 
Euclidean (perpendicular) distance from v(x) to £(v(x),x). For r G [l,oo) let £ r (v(x),x) be the line parallel to 
e(x) such that 

d(v(x),£ r (v(x),x)) = rd(v(x), £{v{x), x)) and d{£{v{x),x),£ r {v(x),x)) < d(v(x), £ r (v(x), x)). 

Let T r (x) be the triangle similar to and with the same orientation as T{y 3 ) having v(x) as a vertex and £ r (v(x), x) 
as the opposite edge. Then the proportional-edge proximity region N PE {x) is defined to be T r (x) n T(y 3 ). 

Furthermore, let £i(x) be the line such that &(x) n T(y 3 ) ^ and r d(y i: &(#)) = d(yi,£(yi, x)) for i = 1,2, 3. 
Then ri(x) n R(yt) ={z£ i?(y t ) : z)) > d(yi,^(x)}, for i = 1, 2,3. Hence rj(x) = Uti( r U^) n i2(y0). 

Notice that r > 1 implies a; G N PE (x) and x G r^x). Furthermore, lim r _ >00 N PE (x) = T(J^ 3 ) for all x G 
r(y 3 ) \ ^3, and so we define N™(x) = T(y 3 ) for all such x. For x G 3^ 3 , we define N r PE {x) = {x} for all 
r G [l,oo]. Then, for x G R(yi), lim I ._ >00 T J (x) = T(^ 3 ) \ {yj,yfc} for distinct and fc. 

Notice that Xi ? ~ P, with the additional assumption that the non-degenerate two-dimensional probability 
density function / exists with support in T(y 3 ), implies that the special cases in the construction of N PE — X 
falls on the boundary of two vertex regions, or at the center of mass, or X G 3^3 — occur with probability zero. 
Note that for such an F, N PE (x) is a triangle a.s. and r^(x) is a convex or nonconvex polygon. 



Figure 1: Construction of proportional-edge proximity region, N PE (x) (shaded region) for an x £ -R(yi). 




Figure 2: Construction of the Ti-region, 2 (x) (shaded region) for an x £ i?(yi). 

3.6 Relative Edge Density of the Underlying Graphs of Proportional-Edge PCDs 

Let X n = {Xi, X2, ■ ■ ■ , X n } be a sample from a distribution F with support in T(J%). Let D n (r) be the 
proportional edge PCD with vertex set V = X n and arc set A defined by (Xi,Xj) £ A Xj £ N PE (Xi). 

Consider the underlying graphs of the data-random PCD, D n (r). Recall that XiXj £ £ an d iff Xj £ N PE (Xi) n 
r x (X t ,N r PE ) and X^ £ £ OI iff Xj £ N r PE (X t ) U T 1 (X t ,N r PE ). 

Let hff(r) := I(XiXj £ £ a „d) = I(Xj £ N PE (X,) n If (X,)) and hg(r) := I(X 4 X, e for) = € 
N PE (Xi) U (^)) for i ^ j. The relative edge density pf^ d (r) := p an d(D n (r)) depends on n explicitly, and on 
F and N PE implicitly. The expectation E \p* (r)l , however, is independent of n and depends on only F and 
N PE . Let Pand (F,r) := E [hff(r)] and iw(F,r) := Cov [h% d (r),hlf(r)}. Then 

0<E[p- d (r)] -E[/^ d (r)] < 1. (7) 

The variance Var [p^ nd (r)] simplifies to 

< Var [pT d (r)] = -T^Var [^ ld (r)] + ^— ^Cov [^ 2 ld (r), ^ d (r)] < 1/4. (8) 
By Theorem it follows that 

( P T d (r) - jw(F, r)) A (0, 4 ^ and (F, r)) (9) 



provided v an d(F, r) > 0. The asymptotic variance of /c™ d (r) is 4 v an <i(F, r) and depends on only F and N PE . 
Thus we need determine only p an d(F, r) and v aTL A{F, r) in order to obtain the normal approximation 

and/ s approx . , / , . 4f and (F, r)\ 
Pn W ~ ■N [Pa.nd(F, r), . (10) 



The above paragraph holds for p° r (r) = p r(-DnW) a l so with p^ d (r) is replaced by p™(r), h^ 2 d (r) and hf 3 ld (r) 
are replaced by h° 2 (r) and h° 3 (r), respectively. 

For r = 1, N P = E 1 (x) n r\ =1 (x) = £(v(x),x) which has zero R 2 -Lebesgue measure. Then we have 

E [ P T d (r = 1)] = E [h\f{r = 1)] = Mand (r = 1) = P(X 2 G A^(Xi) H IT^i)) = 0. 

Similarly P({X 2 ,X 3 } C Np^-(Xi) n r^ =1 (Xi)) = 0. Thus, !/ and (r = 1) = 0. Furthermore, for r = oo, 
N P = E ™{x) n ri =0 °(a;) = T(y 3 ) for all i G T(y 3 ) \ y 3 . Then 

E [p™ d (r = oo)] = E [^ d (r = oo)] = /x and (r = oo) = P(X 2 G JVj.S°°(Xi) n rr°°(*i) = P(A 2 G TQ> 3 )) = 1. 

Similarly, P({A 2 , A 3 } C A^f 30 ^) n r^pfi)) = 1- Hence v and {r = oo) = 0. Therefore, the CLT result in 
Equation (|10[) does not hold for r G {l,oo}. Furthermore, p^ nd (r = 1) = a.s. and p„ (r = oo) = 1 a.s. For 
r G (l,oo), since /if 2 ld (r) = I(X 2 G N PE (Xi) n rj(.X'i) is the number of edges in the AND-underlying graph, 
/if 2 d (r) tends to be high if the intersection region is large. In such a case, hi 3 d (r) tends to be high also. That is, 
hf 2 ld (r) and hf 3 (r) tend to be high and low together. So, for r G (1, oo), we have i^nd^, r ) > 0- See a l so Figure 
[3] (right) and Appendix 1 . 

Forr = l,Np" E 1 (x)UT r 1 = 1 (x) has positive R 2 -Lebesgue measure. Then P({X 2 , X 3 } C N p ^ 1 (X 1 )UTl =1 (X 1 )) > 
0. Thus, i/ or (r = 1) ^ 0. On the other hand, for r = oo, A™(Ai) U r£=°°(Xi)) = T(y 3 ) for all X 1 G T(^ 3 ). 
Then 

E [p° r (r = oo)] = E [hH(r = oo)] = P(X 2 G iVS 00 ^) U rr°°(^i)) - Mor(r = oo) = P(X 2 G = 1. 

Similarly, F({X 2 ,X 3 } C N P %°°(Xi) U rj=°°(Xi)) = 1. Hence i/ or (r = oo) = 0. Therefore, the CLT result 
for the OR-underlying case does not hold for r = oo. Moreover p° r (r = oo) = 1 a.s. For r € [l,oo), since 
h° T 2 (r) = I(X 2 G N PE (Xi) U V[(Xi) is the number of edges in the OR-underlying graph, h° 2 (r) tends to be high 
if the union region is large. In such a case, h° 3 (r) tends to be high also. That is, h° 2 (r) and h° 3 (r) tend to be 
high and low together. So, for r G [1, oo), we have v m (F, r) > 0. See also Figure [3] (right) and Appendix 2. 

Remark 3.5. Relative Arc Density of Proportional-Edge PCDs: 

Let hij(r) := I((Xi,Xj) G A) = l(Xj G N PE {Xi)) for i ^ j and the relative arc density p n (r) := p{D n {r)). 
Let p{r) :=E[p„(r)] and4^(r) := Cov [hi 2 {r), hi 3 (r)]. By Theorem 12. 2i we have 

VE( Pn (r)-p(r))-^Af(0,4v(r)) (11) 
pr ovided v(r) > 0. The explicit forms of asymptotic mean p(r) and variance 4 v(r) for uniform data are provided 



Cevhan et al.l (|2006l) . □ 



4 Asymptotic Distribution of Relative Edge Density for Uniform 
Data 

Let Xi *~ U(T(~y 3 )) for i = 1, 2, . . . , n, where U(T(y 3 )) is the the uniform distribution on the triangle T(y 3 ). 

We first present a "geometry invariance" result which will simplify our subsequent analysis by allowing 
us to consider the special case of the equilateral triangle. Let p™ d {r) := p^d{U(T{y 3 )), r) and p° r (r) := 
PoI (U(T(y 3 )),r). 

Theorem 4.1. Geometry Invariance: Let y 3 = {yi,y 2 ,y3} C R 2 be three non-collinear points. For i = 

1,2, ... ,n, let Xi ~ U(T(y 3 )). Then for any r G [1, oo] the distribution of Pn nd {r) and p°^(r) is independent of 
y 3 , and hence the geometry ofT{y 3 ). 

Proof: A composition of translation, rotation, reflections, and scaling will take any given triangle T a = T(yi, y 2 , y 3 ) 
to the "basic" triangle T b = T((0, 0), (1, 0), (ci, c 2 )) with < Ci < 1/2, c 2 > and (1 - cj) 2 + c\ < 1, preserving 

uniformity. The transformation <p : R 2 —> R 2 given by <p(u, v) = (u+ 1 2 2 c ^ 1 v, v) takes Tb to the equilateral 



triangle T e = T ((0,0), (1,0), (1/2, v3/2)). Investigation of the Jacobian shows that cj> also preserves uniformity. 
Furthermore, the composition of </> with the rigid motion transformations and scaling maps the boundary of the 
original triangle T a to the boundary of the equilateral triangle T e , the median lines of T D to the median lines of 
T e , and lines parallel to the edges of T D to lines parallel to the edges of T e . (A median line in a triangle is the line 
joining a vertex with the center of mass.) Since the joint distribution of any collection of the h^ A (r) and h°j(r) 
involves only probability content of unions and intersections of regions bounded by precisely such lines, and the 
probability content of such regions is preserved since uniformity is preserved, the desired result follows. I 

Based on Theorem 14.11 for our proportional-edge proximity map and the uniform data, we may assume that 
T(ys) is a standard equilateral triangle, T e , with vertices 3^3 = {(0, 0), (1, 0), (1/2, V3/2)}, henceforth. 

In the case of the (proportional-edge proximity map, uniform data) pair, the asymptotic distribution of pf^ d (r) 
and p° r (r) as a function of r can be derived. Recall that p and (r) = E [/if' 2 ld (r)] = P(X 2 G N$, E (Xi) n rj(J?i)) 
and p or (r) = E [hf 2 ] = P(X 2 G N£ E (Xi) U Tl(X{)) are the edge probabilities in the AND- and OR-underlying 
graphs, respectively. 

Theorem 4.2. (Main Result 3) For r e (l,oo), 

(/CV) - p and (r)) j^\v and {r) A Af(0, 1) 

and for r G [1, oo), 

Vn~(Pn r (r)-Por(r)) /V^j A^(0,1). 
where the asymptotic means are 



f (1 - r)(B ^- 1 ^^grrl?) ra " 232 r+128) fOT re [1,4/3), 

101r 5 -801r 4 + 1302r 3 -732 r 2 -536r+672 „ ! A in o /o\ 

216r(r+2)(r+l) IOT 7 fc PV' 3 )' 3 /^! 

r 8 - 13 r 7 +30 r 6 + 148 r 5 -448 r 4 +264 r 3 +288 r 2 -368 r+96 

8r 4 (r+2)(r+l) 
(r 3 +3r 2 -2+2r)(-l+r) 2 
r 4 (r+l) 



for re [3/2,2), 
for r G [2, oo), 



(12) 



Pot(t) 



r 4 (r+l) 

and the asymptotic variances are 



' 47 r 6 -195 r 5 +860 r 4 -846 r 3 -108 r 2 + 720 r-256 f £ M ^/oA 

108 r 2 (r-t-2)(r-f-l) L ' ' 

175 r" -579 , 4 + 1450^ -732 , 2 -536 r+672 for ^[4/33,3) 

216 r (r+2)(r+l) L / * / / * n\ 
3 r 8 -7 r 7 -30 r 6 +84 r 3 -264 r 4 +304 r 3 + 144 r 2 -368 r+96 iv _ c [0/0 OA ^ ' 
8r<t(r+2)(r + l) tOT r e l-V-MJ) 

^tC+V 2 for re [2,oo), 



11 



4=1 



11 



^(^VfWIft). (15) 



i=l 



The explicit forms of $f nd (r) and $° r (?") are provided in Appendix Sections 1 and 2, and the derivations of PandM, 
,y and( T "), p r(?"), and f r( 7 ') are provided in Appendix Section 3. 

The expectation E [/if^M] = Paad( r ) is as in Equation (|12l) ; and E [h° 2 (r)} = Poi (r) is as in Equation (|13l) 
(see Figure 0J Notice that /x an d(r = 1) = and lim r ^ 00 p an d(/") = 1 (at rate 0(r -1 )); and fJ, OI (r = 1) = 37/108 
and limr^oo p or (r) = 1 (at rate 0(r -1 )). 

To illustrate the limiting distribution, for example, r = 2 yields 



V^(pr d (2) - Ma„d(2)) _ / 362880 n ( p ,n A{2) _ 11\ 



v/4^a„d(2) V 58901 V™ V ' 24 



and 



or equivalently, 



0i(p° r (2) - M or(2)) _ / 120960 n ^ or(2) _ 19\ 



^4^ or (2) V 13189 V 24 



- d (2) ap K'° x ('-, J^L) and p « (2) ~ U ( __, _______ 

1 24' 362880 n Hn w V 24' 120960 n 



Figure 3: Result of Theorem 14.21 asymptotic null means p(r), p a nd(f), an d Por{r) (left) and variances &i/(r), 
4i^and(r), and 4^ or (r) (right) for r € [1,5]. Some values of note: /i(l) = 37/216, // an d(l) = 0, and /x or (l) = 
37/108, limr^-oc p(r) = lim r ^ 00 p an d(f) = lim r _ ) . 00 p or (r) = 1, 4y aI1 d(f = 1) = and lim r ^oo 4 i/ and (r) = 0, 
4f r(f = 1) = 1/3240 and lim r _>oo 4 ^ or (r) = 0, and argsup,,^ , 4 z/(r) ps 2.045 with sup re r l oo i 4 z/(r) sa 
.1305, argsup re[l oo] 4i/ and (r) w 2.69 with sup re[l oo] 4 i/ and (r) ps .0537, argsup re [ loo ] 4 ^ or (r) ps 1.765 with 

SU Pre[l : oo] 4 ^orW ~ -0318. 




Figure 4: The edge probabilities p and (r) (left) and p i(r) (right) for r £ [1, 5]. 



By construction of the underlying graphs, there is a natural ordering of the means of relative arc and edge 
densities. 

Lemma 4.3. The means of the relative edge densities and arc density (i.e., the edge and arc probabilities) 
have the following ordering: p an d(r) < p(r) < p r{r) for all r € [1, oo). Furthermore, for r = oo, we have 
Pand(r) — p{r) — p or {r) = 1. 

Proof: Recall that p and (r) = E[^ nd (r)] = P(X 2 e A^(*i) n rf(X x )), p(r) = E[p„(r)] = P(X 2 e 2V£ B (Xi)), 
andp or (r) = E[p° r (r)] = P(X 2 e JV£ B (X0uTJpTi)). And [JV£ B (JTi) PI rj(JTi)] C N r PE {Xi) C [iV£ B (Xi) U rj(JTi)] 
with probability 1 for all r > 1 with equality holding for r = oo only. Then the desired result follows (See also 
FigureE]). ■ 

Note that the above lemma holds for all X^ that has a continuous distribution on T(y s ). There is also a 
stochastic ordering for the relative edge and arc densities as follows. 

Theorem 4.4. For sufficiently small r, p% nd (r) < p n (r) < ST p°^{r) as n — > oo. 

Proof: Above we have proved that p and (?~) < p(r) < p r(»") for all r £ [1, oo). For small r (r < 1.8) the asymptotic 
variances have the same ordering, 4z/ and (r) < Av{r) < Av m {r). Since pf* ld (r), p ra (r), Pn( r ) are asymptotically 
normal, then the desired result follows (See also Figure [3]). I 

We assess the accuracy of the asymptotic normality for finite sample data based on Monte Carlo simulations. 
We generate n X points independently uniformly in the standard equilateral triangle T e . For each data set 
generated, we calculate the relative edge density values for the AND- and OR-undcrlying graphs based on the 



proportional-edge PCD with r = 2. We replicate the above process N mc = 1000 times for each of n = 10, 20, 
and 100. We plot the histograms of the relative edge densities of the AND- and OR-undcrlying graphs using 
the simulated data and the corresponding (asymptotic) normal curves in Figures [5] and H3 respectively. Notice 
that, for r = 2, the normal approximation is accurate even for small n although kurtosis may be indicated for 
n = 10 in the AND-underlying case, and skewness may be indicated for n = 10 in the OR-undcrlying case. 
We also investigate the behavior of the relative edge densities for extreme values of n and r. So we generate 
n = 10 X points and calculate the relative edge densities for r = 1.05 and r = 5. We repeat the above procedure 
N m c = 10000 times and plot the histograms of the relative edge densities in Figures [7] and [SJ which demonstrate 
that severe skewness is obtained for these extreme values of n and r. The finite sample variance and skewness 
may be derived analytically in much the same way as was 4i/ an( j(r) (and Au OT {r)]) for the asymptotic variance. 
In fact, the exact distribution of p^ d {r) (and Pn( r )) is, in principle, available by successively conditioning on the 
values of the Xi. Alas, while the joint distribution of /if| (r), hf^ A (r) (and h™ 2 {r) : h°^(r)) is available (see Figures 
[91 and [P0|) . the joint distribution of {hfj d (r)}i<i<j< n (and {h™(r)}i<i < j< n ), and hence the calculation for the 
exact distribution of p™ d (r) (and p£ r (r)), is extraordinarily tedious and lengthy for even small values of n. 




Figure 5: Depicted are pT d {2) ap -° x J\f for n = 10 ' 20 ' 100 ( left to ri § ht )- Histograms are based 

on 1000 Monte Carlo replicates. Solid lines are the corresponding normal densities. Notice that the axes are 
differently scaled. 




Figure 6: Depicted are p° r (2) ap -° x M (±f, j^§^) for n = 10,20,100 (left to right). Histograms are based 
on 1000 Monte Carlo replicates. Solid lines are the corresponding normal densities. Notice that the axes are 
differently scaled. 

Let 7n(r) be the domination number of the proportional-edge PCD based on X n which is a random sample 
from U(T(y^)). Additionally, let 7^ nd (r) and 7™(r) be the domination numbers of the AND- and OR-undcrlying 
graphs based on the proportional-edge PCD, respectively. Then we have the following stochastic ordering for the 
domination numbers. 

Theorem 4.5. For all r G [l,oo) and n> 1, 7 I ° r (r) < ST 7„(r) < ST 7°™ d (r). 



Figure 7: Depicted are the histograms for 10000 Monte Carlo replicates of pfg d (1.05) (left) and p?g d (5) (right) 
indicating severe small sample skewness for extreme values of r. Notice that the vertical axes arc differently 
scaled. 




Figure 8: Depicted are the histograms for 10000 Monte Carlo replicates of p° T (l) (left) and p°o(5) (right) 
indicating severe small sample skewness for extreme values of r. Notice that the vertical axes arc differently 
scaled. 



Proof: For all x € T(y 3 ), we have [N r PE {x) n T\{x)] C N r PE {x) C [JV^(i)Uf[(i)]. For X ~ U{T(y 3 )), we 
have [N PE (X) n T[ (X)] C N PE (X) C [A^ B (X) U T[ {X)] a.s. Moreover, 7 „(r) = 1 iff X n C N PE (X t ) for some 
i; 7r d ( r ) = 1 iff x n C N PE (X t ) n rj (X,) for some i; and 7° r (r) = 1 iff X n C N PE (Xi) U rj (X,) for some 
i. So it follows that P( 7 £ nd (r) = 1) < P(7„(r) = 1) < P(7° r (r) = 1). Similarly, for all x,y € T(y 3 ), we have 



For 



{[N PE (x) n r{(x)} U [JV£ B (y) n rfty)]) C (JVJ B (x) U iV^(j/)) C ([JVJ B (x) U rj(x)] U [jV^(y) u n(y)] 

x,y # w(r(y 3 )), we have ([iv^(x) n r;(*)] u [jv^(y) n rj(y)]) c (Xj; B (x) u a^(y)) c ([iv^x) u 

FJ(X)] U [N PE (Y) U rj(y)]) a.s. Moreover, 7n (r) < 2 iff C AWX) U iV^(^) for some i ^ j; 7 ^ d (r) < 2 
[X£ £ (X. t ) n rj (X t )] U [N r PE (X 3 ) n r; (X,-)]) for some i ^ j; and 7 ° r (r) < 2 iff X n C ([X£ S (X^) u 



iff x n c 



rj (Xi)] U [Np^Xj) U (X 3 -)] ) for some i j. So it follows that P(7* nd (r) < 2 ) < P(7»M < 2) < P( 7 "(r) < 
2). Since P( 7n (r) < 3) = 1 (fCevhan and Priebej (120051 )). it follows that P(7° r (r) < 3) = 1 also holds since 
P{ln{ r ) < 3) < P(jn(r) < 3). Hence the desired stochastic ordering follows. ■ 

Note the stochastic ordering in the above theorem holds for any continuous distribution F with support being 



in 



T(y 3 ). For r = oo, we have 7° r (r) = 7 „(r) = 7 „ (r) = 1 a.s 



5 Multiple Triangle Case 

Suppose y m is a finite set of m > 3 points in M 2 . Consider the Dclaunay triangulation (assumed to exist) of y m . 
Let Ti denote the i th Delaunay triangle, J m denote the number of triangles, and Cr (3^m) denote the convex hull 

of y m . For Xi ~ U(CH(y m )), i = lj 2, . . . ,n, we construct the proportional-edge PCD, D n ^ m (r), using N PE {-) 




Figure 9: The plots for the joint distribution of hff(r), hff(r) for r £ [1, 5]. Plotted are P(h^ d (r), hff(r) = 
(0,0)) (left), P(htf(r),hff(r) = (1,0)) = P(hff(r),hff(r) = (0,1)) (middle), and P(h% d (r), hff(r) = (1,1)) 
(right). 





Figure 10: The plots for the joint distribution of h° 12 {r), hf 3 (r) for r G [1,5]. Plotted are P(/i° r 2 (r), hf 3 (r) = (0,0)) 
(left), P(h<&(r),V&(r) = (1,0)) = P(h? 2 (r),h%(r) = (0,1)) (middle), and P(h? 2 (r),h%(r) = (1,1)) (right). 



as described in Section [3. 5 1 where for Xi E Tj, the three points in y m defining the Delaunay triangle Tj are used 
as y^] . We investigate the relative edge densities of the underlying graphs based on the proportional-edge PCD. 
We consider various versions of the relative edge density in the multiple triangle case. 

5.1 First Version of Relative Edge Density in the Multiple Triangle Case 

For J m > 1, as in SectionEH let pj nd (r) = 2 |£ and | /(n (n - 1)) and p°*{r) = 2 |£ or | /(n (n - 1)). Let £^ d be the 
number of edges and p™ d (r) be the relative edge density for triangle i in the AND- underlying case, and £°q and 
p°* (r) be similarly defined for OR-underlying case. Let be the number of X points in Tj for i = 1, 2, . . . , J rn . 
Letting Wi = A(Ti)/A(CH(y m )) with A(-) being the area functional, we obtain the following as a corollary to 
Theorem 



(16) 



Corollary 5.1. For r € (l,oo), the asymptotic distribution for pf „ (r) conditional on y m is given by 

(p?#(r) - Pand(m, r)) A N (0, 4Zw(m, r)) , 
as n oo, where p a nd{m, r) = p an d{r) (j2i=i w f) and 



Vand{m, r) 



Vand{r) X^^H + (P<™d(r)) 2 Yl' 



i=l 




with p an d(r) and is an d(r) being as in Equations (|12|) and (|14[) , respectively. The asymptotic distribution of pj r n (r) 
with r £ [l,oo) is similar. 

The proof is provided in Appendix 4. By an appropriate application of the Jensen's inequality, we see 
that J2i=i w i — \ J2i=i w i) ■ S° the covariance above is zero iff is an d(r) = and J2i=i w i = (j2i=i w i J > so 
asymptotic normality may hold even though v an d(r) = 0. That is, p1 nd {r) has the asymptotic normality for 



r G {1, oo} also provided that J2i=i w i > 



The same holds for the OR-underlying case (for r = oo). 



5.2 Other Versions of Relative Edge Density in the Multiple Triangle Case 

Let ~rV) := J2 Zt~l) P " d{r) - Tlim S " nd(r) = P?n " (r) ' SinC<3 S " ld(r) = £ Zt-l) ^ = 
V Jm 2 |£ and | o I r I 

, " = 4^4 - ^ d M- Similar1 ^ we have s « w = >w 

n (n — 1 ) n (n — 1) 

Furthermore, let S^ nd :— Y^,i=i w 1 Pp" d ( r ) where Wi is as in Section l5~Tl So Zf* ld a mixture of the p?? d (r) values. 
Since the p^ d (r) are asymptotically independent, S and (r), pf\ d (r) are asymptotically normal; i.e., for large n their 
distribution is approximately J\f (p an d( m , r ), 4 ? a nd( m : r )/ n )- A similar result holds for the OR- underlying case. 

In Section I5~TI the denominator of pj nd (r) has n(n — l)/2 as the maximum number of edges possible. How- 
ever, by definition, given the rii we can at most have a graph with J m complete components, each with order 
rij for i = 1,2,..., J m . Then the maximum number of edges possible is n t := Yli=i n i ( n i ~ which sug- 

\£ I y^ J _l" l£ a " d l 

gests another version of relative edge density, namely, p^ d n {r) := Then pf? d „(r) 



n t n t 

V Ui ^ - ^ pf nd (r). Since Wi t"*" 1 ) > for each i, and V ( "' ~ 1) = 1, ^ (r) is a mixture of the p a " d (r). 
^ — ' 2n t w ' 2 n t ' w 

Then E [p/j^r)] = p a nd{r), and the asymptotic variance of Pxi d n (r) is 



2\ 




4. and (r) / E u '< + 4 ^nd(r)) 2 £ 




1 

72 

A similar result holds for the OR-undcrlying case also. 

Theorem 5.2. The asymptotic distribution for pff d n (r) conditional on y m for r G (l,oo) is given by 

£ 



fo{pu%( r )-Pa,nd{m,r)) -^Af(OAvand(m,r)) , (17) 



as n -> oo, where p a nd{m,r) = p an d(r) and i> a nd(m,r) 



id 



2 

2 ' 



with Pandir) ana 



Vand{f) being as in Equations (|12[) and (|14p . respectively. The asymptotic distribution of p°j (r) r G [l,oo) 
«s similar. 

The proof is provided in Appendix 5. Notice that the covariance t'and^, ?*) is zero iff is an d(r) = 0. The 
OR-underlying case is similar. 

Remark 5.3. Comparison of Versions of Relative Edge Density in the Multiple Triangle Case: Among 
the versions of the relative edge density we considered, S a ! nd (r) = pj nd (r) for all n > 1, and S and and Pj nd (r) are 
asymptotically equivalent (i.e., they have the same asymptotic distribution). However, pj nd (r) and pf} d n {r) do not 
have the same distribution for finite or infinite n. But we have Pj nd (r) = n ^l-^ pf} d n (r) and since w 1 < 1> ^ 

follows that p and (m, r) < p an d(m, r) = p an d(r). Furthermore, since n ^li) = Yn=i ^("'-"l^ — > Yh=\ w i as n * ~* 
oo, we have linin^oo Var[y/np^(r)} = (j2i=i w f) linV->oo Var[^/np^(r)]. Hence v and (m, r) > V and (m,r). 
Therefore, we recommend Pj nd (r) for use in spatial pattern analysis in the multiple triangle case. Moreover, 
asymptotic normality might hold for Pj nd (r) even if v and {r) — 0. □ 



5.3 Extension to Higher Dimensions 

The extension to M. d for d > 2 is straightforward. Let y>d+i = {yi,Y2, ■ • ■ ,Yd+x} be d + 1 non-coplanar points. 
Denote the simplex formed by these d + 1 points as 6(3^z+i)- A simplex is the simplest polytope in Mr having 
d+ 1 vertices, d (d + l)/2 edges and d + 1 faces of dimension (d — 1). For r G [1, oo], define the proportional-edge 
proximity map as follows. Given a point x in ©(34j+i)i let y := argmin ye ^ tJ+1 volume(Q y (x)) where Q y (a;) is 
the polytope with vertices being the d(d+ l)/2 midpoints of the edges, the vertex y and x. That is, the vertex 
region for vertex v is the polytope with vertices given by v and the midpoints of the edges. Let v(x) be the 
vertex in whose region x falls. If x falls on the boundary of two vertex regions or at the center of mass, we assign 
v(x) arbitrarily. Let ip(x) be the face opposite to vertex v(x), and rj(v(x),x) be the hyperplane parallel to f{x) 



which contains x. Let d(v(x),r](v(x), x)) be the Euclidean distance from v(x) to r](v(x), x). For r £ [l,oo), let 
r] r (v(x),x) be the hypcrplanc parallel to ip(x) such that 

d(v(x),r] r (v(x),x)) = rd(v(x),T)(v(x),x)) and d(r](v(x),x),r] r (v(x),x)) < d(v(x), r/ r (v(x), x)). 

Let & r (x) be the polytope similar to and with the same orientation as © having v(x) as a vertex and rj r (v(x), x) 
as the opposite face. Then the proportional-edge proximity region N PE (x) :— & r (x) n ©(J^j+i)- Furthermore, 
let Q(x) be the hyperplane such that Q(x) n ©(3^+i) 7^ and r d(yi, d(x)) = d(yi,r](yi,x)) for i = 1,2, . . . , d+ 1. 
Then rfta;)nii(yO = {z e ii(y<) : d( yil *?(y<, «)) > <*(*•,&(*)}, for » = 1,2,3. Hence T{(x) = \jj±l{TUx)nR(yi)). 
Notice that r > 1 implies .t € Np E (x) and x <E IT (a:). 

Theorem 14 . 1 1 generalizes . so that any simplex & in M. d can be transformed into a regular polytope (with edges 
being equal in length and faces being equal in volume) preserving uniformity. Dclaunay triangulation becomes 
Delaunay tessellation in R d , provided no more than d + 1 points being cospherical (lying on the boundary of 
the same sphere). In particular, with d = 3, the general simplex is a tetrahedron (4 vertices, 4 triangular faces 
and 6 edges), which can be mapped into a regular tetrahedron (4 faces are equilateral triangles) with vertices 
(0, 0, 0) (1, 0, 0) (1/2, VS/2, 0), (1/2, V3/4, VS/2). 

Asymptotic normality of the [/-statistic holds for d > 2 in both underlying cases. 



6 Discussion and Conclusions 



In this article, we demonstrate that the relative edge density of random graphs and relative arc density of random 
digraphs are one-sample [/-statistics of degree 2. Then, we specify the conditions under which the asymptotic 
normality of the relative densities follows for the random graphs and digraphs. We consider the asymptotic 
distribution of the relative edge density of the underlying graphs based on (parametrized) proportional-edge 
proximity catch digraphs (PCDs). In particular, we consider the AND- and OR- underlying graphs based on the 
proportional-edge PCD; and derive the asymptotic distribution of the relative edge density using the central limit 
theory of [/-statistics. We compute the asymptotic mean and variance of the limiting normal distribution for 
uniform data based on detailed geometric calculations. Moreover, we compare the asymptotic distributions of the 

relative edge densities of the underlying graphs and of the relative arc density of the PCDs. 

The PCDs have applications in classification and spatial pattern analysis. ICevhan et al.1 (|2006l ) showed that 
the relative (arc) density of the proportional-edge PCDs is a [/-statistic and used it for testing bivariate spatial 
patterns. The relative edge densities of the underlying graphs based on this PCD can be employed for the same 
purpose. More specifically, the relative edge densities can be employed for testing the complete spatial randomness 
(CSR) of two or more classes of points against the segregation or association of the points from the classes. CSR 
is roughly defined as the lack of spatial interaction between the points in a given study area. In particular, the 
null hypothesis can be assumed to be CSR of X points, i.e., the uniformness of X points in the convex hull of y 
points. Segregation is the pattern in which points of one class tend to cluster together, i.e., form one-class clumps. 
On the other hand, association is the pattern in which the points of one class tend to occur more frequently 
around points from the other class. Under the segregation alternative, the X points will tend to be further away 
from y points and under the association alternative X points will tend to cluster around the y points. Such 
patterns can be detected by the test statistics based on the relative edge densities, since under segregation we 
expect them to be smaller, and under association they tend to be larg er. Moreover, the underlying graphs can 
also be used in pattern classification as outlined in lPriebe et al.l (|2003al ). Finally, the methodology described here 
is also applicable to PCDs in higher dimensions. 
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APPENDIX 

Appendix 1: The Asymptotic Variance of Relative Edge Density for the AND- 
Underlying Graph Version: 

The variance of hf^ir), denoted as Var an( j(?') = Var [/ii^M] ■ 1S as follows: 

Var and (r) = <\ d (r)I(r e [1,4/3)) + < 2 d (r)I(r 6 [4/3,3/2)) + < 3 d (r)I(r e [3/2,2)) +^/(r)I(r e [2, M )) 

V, and/ \ _ (5 r 6 -153 r 5 +393 r 4 -423 r 3 -54 r 2 +360 r-128)(447 r 4 -261 r 3 +54 r 2 +5 r 6 -153 r 5 +360 r-128) 

wnere Lp ll (r) - - 2916r 4 (r+2) 3 (r+l) 2 > 

and/ ^ _ (101r 5 -801r 4 + 1302r 3 -732r 2 -536r+672)(1518r 3 -84r 2 -104r+101r 5 -801r 4 +672) 
</>l,2 l r J — 46656r 2 (r+2) 2 (r+l) 2 ' 

and/ \ _ (r 8 -13r T +30 r 6 + 148 r 5 -448 r 4 +264 r 3 +288 r 2 -368 r+96) (22 r 6 + 124 r 5 -464 r 4 +r 8 - 13 r T +264 r 3 +288 r 2 -368 r+96) 
^1,3 l r J — ~~ 64rS(r+2) 2 (r+l) 2 ' 

^iV) = (r5+r4 " 3r3 ~ 3 r2 r ty + " 1 2 ) l (3r3+3r2 ' 6 r+2) ■ Note that Var and (r = 1) = and lim™ Var and (r) = (at 
rate 0(r~ 2 )), and argsup re r loo ) Var an d(r) w 2.1126 with sup Var an d(r) = 0.25. 
The asymptotic variance for the AND-undcrlying version is 

n 

^ and (r) := Cov [/^ d (r), htf(r)] = £ Ifr) 

i=i 

where 

^(r) = - 58320 (2f . 2 + 1)( t + 2)2(r + 1) 3 r e ((r " D^r 19 + 8748 r" + 44456r" + 140328 r" + 121371 r^ 

- 412117 r 14 - 27145 r 13 - 4503501 r 12 + 1336147 r 11 + 10640999 r 10 - 982009 r 9 - 6677105 r 8 - 2274458 r 7 

- 1150162 r 6 + 249126 r 5 + 1232530 r* + 1234372 r 3 + 226776 r 2 - 184944 r - 81920)) 

$T d M = — — — - (486 r 21 + 3402 r 20 - 269 r 19 - 45155 r 18 - 118850 r 17 + 443518 r 16 

2 v; 116640 (2r 2 + l)(r + 2) 2 (r + l) 3 r 6 v 

+ 3251855 r 15 - 13836295 r 14 + 13434672 r 13 + 11140788 r 12 - 27667544 r 11 + 13293088 r 10 + 7159710 r 9 - 

13013598 r 8 + 4185440 r + 3262952 r 6 + 586636 r 5 - 1616444 r 4 - 680120 r 3 - 55952 r 2 + 219936 r + 49152) 

$T d (r) = ; — 5 rrr, rrr-jr (486 r 21 + 3402 r 20 - 269 r 19 - 45155 r 18 - 118850 r 17 + 443518 r 16 

3 w 116640 (2r 2 + l)(r + 2) 2 (r + l) 3 r 6 v 

+ 2751855 r 15 - 13736295 r 14 + 18084672 r 13 + 8770788 r 12 - 43009544 r 11 + 24604048 r 10 + 27137438 r 9 - 30889822 r 8 
- 2832544 r + 11101160 r 6 - 4168820 r 5 + 2364868 r 4 + 2305864 r 3 - 3041936 r 2 + 219936 r + 49152) 



*r d M = - 58320{r + mr2 J ){2r2 + 1){r + 1 ys rG 0632 , 22 + 25632 , 24 - 60328 - 441888 r 19 + 1353430 r 18 

- 297666 r 17 - 4791125 r 16 + 12849927 r 15 - 10894618 r 14 - 26295324 r 13 + 62283823 r 12 - 2280753 r 11 - 81700012 r 10 
+32551926 r 9 +39974410 r 8 -11284026 r 7 -5806580 r 6 -9167580 r 5 -2004944 r 4 +4646688 r 3 +1931776 r 2 -489024 r-98304) 

■&T d (r) = C*V) = - 1 — — — (3632r 22 +25632r 21 -49432r 20 - 36 4992r 19 +958940r 18 



58320 (r + 2) 3 (2r 2 + l)(r 2 + l)(r + l) 3 r 6 
.8 r 16 + 5424126 r 15 - 23566328 r 14 + 238: 
- 8239197 r 9 - 30178496 r 8 + 27901506 r - 4936170 r 6 + 61038 r 5 + 4719720 r 4 - 5513952 r 3 + 340736 r 2 + 23328 r + 65536) 



1167012 r 17 + 1200518 r 16 + 5424126 r 15 - 23566328 r 14 + 23837088 r 13 + 11797395 r 12 - 41623065 r 11 + 39261953 r 10 



i9? nd W = 77777777 77777 777 777"^ (1562 r 21 - 11142 r 20 - 103099 r 19 + 2105697 r 18 - 9774118 r 17 + 

w 466560 (r + 2) 3 (2 r 2 + l)(r 2 + l)(r + l) 3 r 5 v 

10220280 r 16 + 27825711 r 1S - 69243129 r 14 + 81624200 r 13 - 76052574 r 12 - 65530400 r 11 + 262451196 r 10 - 178092280 r 

- 69106464 r 8 + 158439568 r 7 - 97568688 r 6 + 12246288 r 5 + 17591952 r 4 - 21111616 r 3 + 15628032 r 2 - 2545664 r + 993024) 



$8 nd (r) = -7777777 7777-7 7777-7 777 777-77 (2r 26 -30r 25 -2395r 23 +281r 24 +8770 r 22 +29 5 28 r 21 - 

8 w 1920(r + 2) 3 (r 2 + l)(2r 2 + l)(r+l) 3 r 10 v 

245667 r 19 + 2066216 r 18 - 5313494 r 17 - 1589216 r 16 + 18512684 r 15 - 18946136 r 14 - 2665248 r 13 + 22789! 
32987760 r 11 + 20482512 r 10 + 13109584 r 9 - 28084416 r 8 + 17326976 r - 3864576 r 6 - 4579328 r 5 + 6666240 r 4 - 

. O . . 



-268053 r 20 + 

+ 22789584 r 12 - 
4579328 r 5 + 6666240 r 4 - 3576320 r 3 
+ 635904 r 2 - 116736 r + 61440) 



^" d(r) = " 1920 (r + 2) 3 (, 2 + l)(2r- 2 + l)( r + l) 3 r^° (2 ^-30 r 25 -2395 r 23 281 r 24 + 8258 r 22 + 3 1064 r 21 -262677 r 20 + 
225443 r 19 + 2052136 r 18 - 5219030 r 17 - 1608928 r 16 + 18337836 r 15 - 18837080 r 14 - 2598688 r 13 + 22736336 r 12 - 
32858736 r 11 + 20384720 r 10 + 12930896 r 9 - 27988416 r 8 + 17416832 r 7 - 3862784 r 6 - 4575488 r 5 + 6638848 r 4 - 3603200 r 3 

+ 640512 r 2 - 107520 r + 63488) 

*iSV) = ' 1920 (r + 2 ) 3 (,-lK r + l ) 3 (2r 2 -lK0 (2r 25 +307 r 23 -32r 24 -2612r 22 + 11572r 21 +21934r 20 - 32 8 8 67r 19 + 

524994 r 18 + 2446870 r 17 - 8676180 r 16 - 437020 r 15 + 36944680 r 14 - 40677696 r 13 - 44860384 r 12 + 106256352 r 11 - 
15515040 r 10 - 98636848 r 9 + 66358080 r 8 + 27142272 r - 42614272 r 6 + 7781120 r 5 + 7327232 / - 3388672 r 3 + 

430592 r 2 - 171008 r + 63488) 



i?n W = ; — 5 77 77-77 (30r + 90r - 127 r - 621 r + 320r 9 + 1568r 8 - 858 r - 1370 r + 909 r & + 

15 (2H — l)(r + iyr w 

295 r 4 - 292 r 3 + 44 r 2 + 6 r - 2) 

andlj = [1,2/VS), X 2 = [2/V5.6/5), Zj, = [6/5,\/5-l), I 4 = [V5 - 1, (6 + 2 y/2)/7), X 5 = [(6 + 
2 V2)/7, 4/3), Z 6 = [4/3, (6 + v/l5)/7), X 7 = [(6 + \/l5)/7, 3/2), I 8 = [3/2, (1 + v/5)/2), Zg = [(1 + V5)/2, 1 + 
1/V2), I 10 = [1 + 1/V2,2), In = [2,oo). See Figured Note that Cov and (r = 1) = and lim^oo i/ and (r) = 
(at rate 0(r~ 2 )), and argsup r6 j li00 ) i>and(>") ~ 2.69 with supz^andM w .0537. 

Appendix 2: The Asymptotic Variance of Relative Edge Density for the OR-Underlyin^ 
Graph Version: 

The variance of h^if), denoted as Var or (r) = Var [h° r 2 (r)], is as follows: 

Var or (r) = ^ (r)I(r G [1,4/3)) + ^ 2 (r)I(r G [4/3,3/2)) + ^(r)I(r G [3/2, 2)) + ^ 4 (r)I(r G [2, oo)) 

, or / \ _ (47r 6 -195r 5 +860r 4 -846r 3 -108r 2 +720r-256)(752r 4 -1170r 3 -324r 2 +47r 6 -195r 5 +720r-256) 

wnere y x>1 [r) - iT!M7^Tr+2>^7+T? ' 

or I \ _ (175 r 5 -579 r 4 + 1450 r 3 -732 r 2 -536 r+672) (1234 r 3 -1380 r 2 -968 r+175 r 5 -579 r 4 +672) 
Pl^V) — 46656 r 2 (r+2) 2 (r+l) 2 ' 




Figure 12: i> and (r) = Cov [hff(r), h^ d (r)] (left) and v OT (r) = Cov [hf 2 (r), hf 3 (r)} (right) as a function of r for 
re [1,5]. 



or / \ (3 r 8 -7 r 7 -30 r 6 +84 r 5 -264 r 4 +304 r 3 + 144 r 2 - 



vTAir) = 2 (r +r ; s \;^ r ~ 1) - See FigurcUU 

Note that Var or (r = 1) = 2627/11664 and lim 
1.44 with supVar or (r) « .25. 



368 r+96) (-22 r 6 + 108 r 5 -248 r 4 +3 r 8 -7 r 7 +304 rf+1 
64r 8 (r+2) 2 (r+l) 2 



44r 2 -368r+96) 



■r— > oo 



Var or (r) = (at rate 0(r 4 )), and argsup rg [ l oo - ) Var 

~ supVar or (r) pa .25. 
The asymptotic variance for the OR-undcrlying 



lg version is 



v OT (r) := C 



where 

$?(r) = - 

5576151 r 16 + 2978697 r 15 - 



58320 (r 



ll 

'■ov[hH(r),hT 3 (r)]=J2^T(r)K 
i=i 

1W „ 9 \ w rr^-fr (1458r 22 + 13122r 21 +50731 

+ l)(2r 2 + l)(r + l) J (r + 2) 3 r 6 v 

7 r 15 - 33432692 r 14 + 37427862 r 13 + 15883834 r 12 - 60944766 r 11 + 49876417 r 10 - 1 
7 - 10290256 r' - 2234754 r 5 + 7085471 r 4 - 5608569 r 3 + 1645826 r 2 - 132876 



1754523 r 9 - 
r + 30824) 



84225 r 19 -19193 r 18 -1823223 r 17 + 

r" + 37427862 r 13 
36606859 r 8 + 32338215 r - 10290256 r 6 - 2234754 r 5 

*?(r) = *f(r) = - 116640 (r2 + 1)(2r2 ; i)(r + 1) 3 (r + 2) 3 r a (1458 r^+13122 ^+62825 r--175011 r 19 + 156014 , 18 - 
3300900 r 17 + 11053023 r 16 + 5055135 r 15 - 67685050 r 14 + 75243552 r 13 + 33155180 r 12 - 120628524 r 11 + 99831906 r 10 - 
4883958 r 9 -74801558 r 8 +64360782 r 7 -19812000 r 6 -3667716 r 5 +14541630 r 4 -11254002 r 3 +3070468 r 2 -413208 r+28880) 

#?(r) = - r— t^t^ rr-. rr^ -t^ (972 r 24 + 8748 r 23 + 29 5 90 r 22 - 149 106 r 21 - 36820 r 20 - 

58320 (r 2 + l)(2r 2 + l)(r 2 - 2)(r + 2) 3 (r + l) A r b v 

986280 r 19 +5942884 r 18 +2883672 r 17 -47189711 r 16 +43450125 r 15 +85975304 r 14 -156173934 r 13 +27378901 r 12 + 123606417 r 

-152209261 r 10 +64653597 r 9 +56621894 r 8 -88962768 r 7 +43754559 r 6 -5940597 r 5 -13006396 r 4 +17019366 r 3 -7037340 r 2 + 

413208 r- 28880) 

C(r) = ' 58320(r 2 + l)(2r 2 + l)(r + l)3( r + 2 ) 3 r 6 (972 r 22 +8748 r 21 
3745643 r 16 + 4573731 r 15 - 29416804 r 14 + 26163354 r 13 + 19600850 r 1 ' - 

"' + 2292994 r 2 + 14580 r - 1944) 



li 



22237963 r 8 



+31534 r 20 -131610 r 19 +261546 r 18 -1552026 r 17 + 

7 31 r 15 - 29416804 r 14 + 26163354 r 13 + 19600850 r 12 - 43126062 
+ 26778663 r 7 - 9107024 r 6 - 115074 r 5 + 3136927 r 4 - 5055609 i 



„ r 11 + 31497249 r 10 - 7381467 r 9 - 
115074 r 5 + 3136927 r 4 - 5055609 r 3 



C(r) = 233280 (r 2 + l)(2r 2 + l)(r + l)3( r + 2)^ r 22 -7290 r 21 ^ 
1844321 r 16 - 33347697 r 15 + 80028903 r 14 - 29292735 r 13 - 98093906 r 12 



181459 r 20 + 102 4 401 r 19 -2691213 r 18 +3921057 r 17 
r 15 + 80028903 r 14 - 29292735 r 13 - 98093906 r 12 + 125034492 r 11 - 

88057996 r 8 - 



+ 



- 33347697 

- 26383068 r - 12851392 



46658244 r 10 - 57216612 
r 6 + 14179848 r 5 - 8656508 r 4 + 1593828 r 3 + 134136 r 2 - 58320 r + 7776) 



r 9 + 



*™ = 233280 (r + m r> + l)(l> + l)(r + l) H r-l)rS (486^-7776^-174169^ + 1205860^-4656806^ + 
8763566 r 18 +7460036 r 17 -63559490 r 16 +91134324 r 15 + 18516450 r 14 -122708655 r 13 + 18577230 r 12 +80410332 r 11 -19357704 r 10 - 
39129236 r 9 +75311048 r 8 -77449360 r 7 +4053376 r 6 +48283912 r 5 -40690240 r 4 +17736336 r 3 -4315680 r 2 +544320 r-31104) 

*rw = geoFTWTiWwTW (2r24 " 30 r " " 161 r " + 107 r21 + 4137 r2 ° " 10685 r U + 8367 r " 8+ 

78713 r 17 - 450859 r 16 + 697707 r 15 + 517846 r 14 - 3723120 r 13 + 6565124 r 12 - 1468692 r 11 - 8695792 r 10 + 9535720 r 9 - 
6773160 r 8 + 526744 r 7 + 10691376 r 6 - 7797264 r 5 + 1137696 r 4 + 523712 r 3 - 2687872 r 2 + 1701888 r - 245760) 

CO) = , ttt r^T T\ — 77T (2 r 25 - 32 r 24 - 129 r 23 + 236 r 22 + 4 1 57 r 21 - 15610 r 20 + 21289 r 19 + 

960 (2r 2 + l)(r + l) 2 (r + 2) 3 (r 2 + l)r 10 v 

67536 r 18 - 511355 r 17 + 1161830 r 16 - 634128 r 15 - 3001568 r 14 + 9512164 r 13 - 11014136 r 12 + 2344968 r 11 + 7126240 r 10 - 

13850504 r 9 + 14466592 r 8 - 3823216 r - 4018976 r 6 + 5155776 r 5 - 4633984 r 4 + 1959808 r 3 - 244480 r 2 - 3584 r - 1024) 

= 960(2^-1)^ + 2)3^-1)^ + 1)^0 < 2 ^ " 34 ^ " 101 ^ + 433 " + 5400 ^ ~ 26982 ^ + 23049 ^ 
166787 r 17 - 717366 r 16 + 1196092 r 15 + 89468 r 14 - 5130844 r 13 + 12748688 r 12 - 11274744 r 11 - 12243496 r 10 + 
33980568 r 9 - 14886656 r 8 - 19910592 r + 20667776 r 6 - 1262208 r 5 - 5402752 r 4 + 2217088 r 3 - 235776 r 2 - 2560 r - 1024) 

or _ _2_ 180 r 8 - 48 r 7 - 648 r 6 + 396 r 5 + 214 r 4 - 190 r 3 + 39 r 2 - 4 r + 1 
llW ~~ 15 (2r 2 -l)(r + l) 2 r 10 

and = [1,2/VS), X 2 = [2/VS.6/5), I 3 = [6/5,^-1), I 4 = [V5 - 1, (6 + 2 V2)/7), X 5 = [(6 + 
2V2)/7,4/3), I 6 = [4/3, (6 + 7l5)/7), X 7 = [(6 + ^/7,3/2), I 8 = [3/2, (1 + V5)/2), 2g = [(1 + V5)/2, 1 + 
1/V2), 2io = [1 + 1/V2, 2), In = [2, 00). See FigureEl Note that Cov or (r = 1) = 1/3240 and linv^^ u OI (r) = 
(at rate 0(r -6 )), and argsup rG r l oo \ Voiir) ~ 1.765 with sup z/ OT (r) ~ .0318. 

Appendix 3: Derivation of the Asymptotic Mean and Variance for 
Uniform Data 

In the standard equilateral triangle, let yi = (0,0), y 2 = (1,0), y3 = (1/2, \/3/2), Mc be the center of mass, Mi be the 
midpoints of the edges e; for i = 1,2,3. Then M c = (l/2,V3/6), Mi = (3/4, v%/4), Af 2 = (1/4,^/3/4), M 3 = (1/2,0). 
Let X n be a random sample of size n from W(T(^)). For x\ = (u,v), £ r (xi) — r v+r %/3 w— \/3 x. Next, let Ni := £ r (xi)C\e3 
and 7V 2 := n e 2 . 

Appendix 3.1: Derivation of j? an d( r ) an d ^ an d( r ) f° r Uniform Data 
Derivation of yU a nd(f) in Theorem 14.21 

First we find /i a nd(>") for r G (1, oo). Observe that, by symmetry, 

MandW = p(x 2 e Nl E (Xi) n ri(Xi)) = 6P(x 2 g iVj(Xi) n rr(jfi),Xi e r.) 

where T s is the triangle with vertices yi, M3, and Mo- Let £ s (r,x) be the line such that r d(y 1 ,£ s (r,x)) = c((yi,ei), so 
l s {r,x) = \/3(l/r - x). Then if x\ G T s is above then N PE (xi) = T(^ 3 ), otherwise, iV£ B (a;i) C T(^ 3 ). 

To compute ^and(^), we need to consider various cases for N PE (X\) and Ti(Xi) given Xi = (a;, y) G T s . See Figures [T3l 
and !14l For any x — (u, v) G T(y), Ti(x) is a convex or nonconvex polygon. Let £i(r, x) be the line between x and the vertex 
yi parallel to the edge e, such that r d(yi,£i(r, x)) — d(yi, £ r (x)) for i = 1, 2, 3. Then ri(x)niZ(y<) is bounded by £i(r, x) and 
the median lines. For x — (u, v), £i(r, x) = — y/3 x + (v+ Viu) /r, £2(1", x) — (v + Vir (x— 1) + \/3(l — u))/r and £3(7", x) = 

(■\/3(r - 1) + 2 u)/(2 r). For r G [6/5, \/5 - 1), there are six cases regarding T"l(x) and one case for Np E (x). See Figure 
[14] for the prototypes of these six cases of Fi (x,Ny). For the AND-underlying version, we determine the possible types 
of Np E (x{) n ri(a;i) for xi G T s . Depending on the location of xi and the value of the parameter r, Np E (xi) Pi Fi(xi) 
regions are polygons with various vertices. See Figure [T5] for the illustration of these vertices and below for their explicit 
forms. 




Figure 13: The cases for relative position of i s {r, x) with various r values. These are the prototypes for various 
types of N r PE {x\). 



q _ I -J3y+3x fl \ p = / ^3i;-3r+3-33: q\ ^ _ ( Viy -6 r+3- 3 x s/3~(- V3y -3+3 x) \ _ / (V3r+V3-2 y)^ 

y%-3+W3s)\ C - ( (^ r ~^+ 2 v)^ v / 3(3r-3+2 % /3y) \ _ / ^3 y+33; ^(^+33:) ' 

fir J ' ^ 5 — I fir ' fir ' ^ r(i ~ \ fir ' fir 



Pt = (1/2,-^/3/6 (2V3ry + 6rx-3)), and P 2 = (-1/2 + (V3r y + 3 r k)/2, -\/3/6 (-3 + V3r y + 3 r x)) ; 

r = / 1/2 V3(2V3y+6 x-3 r) \ _ / (-2 y^ B -6+6 g+3 r) yg \ _ / V3 y - 3 r+3-3 x V^( 3 r-V3y-3+3 x) \ 

11/) g r 1)21/5 Q r / ' ^ \ ' 6r 1)4 

f 3r-3+2y / 3a V5(3 r-3+2 y^a/) \ _ / r -3+2 y^y V5(3 r-3+2 y^y) \ _ / -^^^3, vg(V3y+3x-3_r) S 

(\/3>2//3 + ra;,0), iV 2 = (v^r y/6 + r x/2, V% (VSy/6 + 3x) r), and 

N 3 = (V3ry/i + 3rx/4,V3(V3y/12 + 3x)r); a nd Qx = ^^+3^-^+3^3+3^ (A'„ + 3r', + ^,-3r + 3-3^ j 



iVi = 



and Q2 



2 v / 3V 2 j/ + 6r 2 :r-3T-+3-2 y^y ^(3 r-3+2 V3y) \ 



Let ^(01,02, • ■ ■ , a,n) denote the polygon with vertices ai, a2, . . . , a„. For r G [l, 4/3) , there are 14 cases to consider 
for calculation of PandfV) in the AND-underlying version. Each of these cases correspond to the regions in Figure 1161 where 
Case 1 corresponds to Ri for i = 1,2, 3, 4, and Case j for j > 1 corresponds to Rj+3 for j = 1,2, . . . , 14. These regions are 
bounded by various combinations of the lines defined below. 

Let £ am (x) be the line joining yi to Mc, then £ a m(x) = y/Sx/3. Let also ri(x) = v3 (2 r + 3 a; — 3) /3, ^(a;) = 
v^/2- y^r/3, ra(aO = (2z-2 + r) v/3/2, r 4 (x) = y/3/2 - y/Sr/4, r 5 (x) = _ 2^^5=11, r6 (a:) = - ^ +3 "' , r 7 (x) = 

fl + r 2 a;-r-a:)v / 3 . . (r 2 i-l + i)v^ , , (r 2 x-l)V3 , . (-2r+2+A)V3 . . (-2 i+r 2 i)v3 
— 1 r i + l I r »\ X ) = ~- T^T ' T *K X ) = ~ r 2 + 2 ' ri °( X ) = ~^ -4+r 2 f11 W = _ ??+2 > 

^12(2;) = — (2a; — r) V3/2, and ri3(a;) = — (—1 + a;) v3/3. Furthermore, to determine the integration limits, we specify 
the a;-coordinate of the boundaries of these regions using Sk for k = 0, 1, . . . , 14. See also Figure [TBI for an illustration of 
these points whose explicit forms are provided below. 

so = 1 - 2r/3, si = 3/2 -r,s 2 = 3/(8 r), s 3 = - 3r + 2 / 2+3 , s 4 = 1 - r/2, s 5 = 2r -£ 2+1 , s 6 = l/(2r), s 7 = 2{ J 2+l) , 

= 9 " 3 6(rt+ r 1 )" 2r , «9 = 1/(^ + 1). sio = , an = 3r/8, s 12 = 8 W = 3/2 - 5r/6, and s 14 = 

r- 1/2 - r 3 /8. 

Below, we compute P(X 2 G N r PE {X 1 ) n r^(Xi),Xi G T 3 ) for each of the 14 cases: Case 1: 

p(x 2 g iVp^Xr) n iuxx), * g r.) = ^ ^ ^ J a^fa = 

(r - 1) (r + 1) (r 2 + 1) 
64 r 6 

where A(^(Gi, JVi, JV 2) G 6 )) = v^/36 (v^y + 3 a:) 2 r 2 - ^(^+ 33: ) 2 



36 r 



y A =(l/2 lV /3/2) 



y A =(l/2,v/3/2) 




Us 



W=ft,0) yi = ((),()) 



Case 2: 



Figure 14: The prototypes of the six cases of T\ (x) for x G T s for r £ [1, 4/3). 



«. A (^(G 1 ,N 1 ,P 2 ,M 3 ,G 6 )) 

L + 4 y J — ^cmw — = 



(9 r 5 + 23 r 4 + 24 r 3 + 24 r 2 + 13 r + 3) (r - 1) 

6 (V 4- I s ! 3 



where A(^(G 1 ,iV lj P 2 ,M 3 ,G 6 )) 
Case 3: 



96 r 6 (r + 1) 

\/3(-4r 3 V3y-12r 3 a; + 2 r 4 y 2 +Ar i \/3y i+6rV+3r 2 +2j 2 +4 \/3</ i + 6s 



P(x 2 e iVp B pfi) n rl(Xi), Xi e r s 



r 7 (i) isg io Jsi2 -/O 



324 r 11 - 1620 r 10 - 618 r 9 + 4626 r 8 + 990 r 7 - 2454 r 6 + 2703 r 5 - 5571 r 4 - 3827 r 3 + 1455 r 2 + 3072 r + 1024 



7776 (r + l) 3 r 6 



where A(&(Gi, G 2 , Qi, Pi, M 3 , G 6 )) = - [V3(-4v / 3r-y - 12 a + 4y 2 + 4r 2 y 2 - 12 r + 9r 2 + 12 r x + 4r 4 y 2 - 12x 2 r 2 
24 r s x + 12 r 4 a; 2 + 8 r 4 V3y x + 12 x 2 + 12 r 2 x + 6 - 8 r 3 ^y + 4 + 4 ^3r 2 y)J / [24 r 2 ] . 



Case 4: 

P(x 2 e n£ b {Xi) n n(Xi),Xi e r«) = 

S5 frjC*) pW ^12 f r 6 (x)\ A{&>(Gl, Ml, L 2 , Ql,P2,M 3 ,G 6 )) , , 

— - — —ay ax = 



r 8 (x) J 8 5 Jr 3 (x) J 8 10 Jr 3 (x) 

512 + 138240 r 7 + 3654 r 12 - 255 r 8 + 43008 r 3 - 12369 r 2 - 86387 r 4 - 193581 r 6 + 148224 r 5 - 100608 r 9 + 94802 r w - 

35328 r 11 ] / [7776 (r 2 + l)V] 

where A{0>(Gi, Mi, L 2 , Qi, P 2 , Ma, G 6 )) = -[V3(6x + 3r 2 - 2 \/3y + 2 x/37- 2 ?/ + 2 r 4 y 2 - 4 r 3 V^y + 4 V^y x + 2r 2 y 2 + 
4 r 4 V3y x-6 x 2 r 2 - 12 r 3 x + 6 r 4 ^ 2 + 6 r 2 x - 3)] / [l2 r 2 ] . 



Case 5: 

p(* 6 w l) n W6 T.)= f r r w + r r w ) ^^^^ ^^ 

Vi S3 •W) isg Jr 6 W / ^(r(^3)) 2 

(177 r 8 - 648 r 7 + 570 r 6 - 360 r 5 + 28 r 4 - 24 r 3 + 174 r 2 + 72 r + 27) (- 12 r + 7 r 2 + 3) 2 

7776 (r 2 + l) 3 r 6 

, »r r n H /, \\ v/3"(-4 r 3 f 3 :c+3 r 2 +6 r 4 v/3i/ z + 9 r 4 x 2 +3 r 4 y 2 + y 2 +2 V3y x + 3 x 2 ) 

where A{&>{Gi, Mi, L 2 , Qi, P 2 , M 3 , G 6 )) = - • 

Case 6: 

P{x 2 e JV^(Xi)nrl(Xi),Xi e r.) = 

r s r 8<3;> ^ A(^(Gi,Mi,Pi,P 2 ,M3,G 6 )) , . 

I L.) J wm 2 dydx= 



€ T.) = 

' s 2 Jr 5 (x) J S3 ^r 2 (a;) 



/ + / / 



\ « o-2 " " °3 " '^l- 1 -; "'2\ J ./ / * \- - / / 

[l37472 r 18 - 952704 r 17 + 2792712 r 16 - 5116608 r 15 + 7057828 r 14 - 7725792 r 13 + 7022682 r 12 - 5484816 r n + 
3631995 r 10 - 2213712 r 9 + 1213271 r 8 - 578976 r + 292518 r 6 - 101952 r 5 + 36612 r 4 - 11664 r 3 + 3051 r 2 - 1296 r + 243] j 

[(I5552r 2 + l) 3 (2r 2 + l) 3 r 6 ] 

U AtOKr* US D D ^ U v^(-4r 3 V3 !/ -12r 3 a: + 3r 2 +6r 4 V3y a: + 9r 4 ;1 ; 2 +3r 4 H 2 +y 2 +2 V^Ji + 3x 2 ) 

where A(^*(Gi, Mi,Pi, P 2 , M 3 , G 6 )j = 1 5^ • 

Case 7: 

P(x 2 e N' M Xi) n r^i), * e r.) = ( f" f W + f° f W ) A(^( Gl ,^ Ql ^,m,, Gb) ) = 

4 (100 r 11 - 408 r 10 + 454 r 9 - 564 r 8 + 283 r 7 - 108 r 6 - 34 r 5 + 204 r 4 - r 3 + 132 r 2 + 26 r + 24) (2 r - l) 2 (r - l) 2 

243 (r 2 + l) 3 r 3 (2r 2 + l) 3 

where A{&>{Gi, Mi, L2, Qi, Pi, Ms, G 6 )) = -[a/3(6x + 3r 2 - 2 %/3y + 2 \/3V 2 y + 2 r 4 y 2 - 4 r^y + 4 x/3y x + 2r 2 y 2 + 

4 r 4 ^3y x-6 xV - 12 r 3 x + 6 r 4 a; 2 + 6 r 2 a; - 3)] j [l2 r 2 J . 
Case 8: 

P(x 2 €ATp B (Xi)nr 1 (Xi),XieT s )= / / +/ / ^ 2 d y dx = 

\Js 12 Jr 6 (x) Js 13 Jr e (x) ) ^K 1 (^3)) 

[(-2 + r) (2369 r 11 - 11342 r 10 + 29934 r 9 - 50340 r 8 + 54056 r - 51824 r 6 + 48320 r 5 - 20864 r 4 - 640 r 3 

- 1280 r 2 + 512 r + 1024)]/ [l5552r 6 J 

where A(5 a (Gi, G 2 , Qi, N s , Ma, Ms, G 6 )) = - [%/3(4^3r 2 i/ - 12 x - 12 r + 5r 2 + 12ri + 4y 2 - 12 xV +4r 2 y 2 + r 4 y 2 + 

2 r 4 V3y x-4 r 3 V3y + 6 - 12 r 3 x + 3 r 4 x 2 + 12 x 2 + 12 r 2 x - 4 \/3V y + 4 \/3y)] / [24 r 2 
Case 9: 

p(x 2 eN^ PB (Xi)nr\(Xi),XieT s )=( f 12 + H H^) A i^Gi 1 M h± 



/raW ; ^(r(^3)) 2 

(49 r 8 - 168 r 7 + 354 r 6 - 528 r 5 + 236 r 4 - 96 r 3 - 224 r 2 + 384 r + 64) (-12 r + 7r 2 + 4 

15552 r 6 



where A{&>{Gi, M x , L 2 , Qi, N 3 , M c , M 3 , G 6 )) = -p3(8^3yz + 4V3r 2 y + 12x + 2r 2 - 12x 2 r 2 - 4r 3 ^3y - 12r 3 x + 

3 r 4 x 2 + r 4 y 2 + 2 r 4 y/?,y x + 12 r 2 x - 6 - 4 \/3y + 4 r 2 y 2 )] / [24 r 2 ] . 
Case 10: 

p(x 2 g iv^pfi) n r;(jCi),jfi g t.) = 

/ r.n rno(-) + p(«) + ri/a rn a («)\ A( ^» (Gi , Mi , L 2 , Qi , JV 3 , L 4 , 1-s , M 3 , G 6 ) ) ^ ^ ^ 

/ A(r(y 3 )) 2 

6144 + 195456 r 6 + 324 r 11 - 76720 r 7 - 801792 r 2 + 217856 r + 946432 r 3 - 239904 r 5 - 275328 r 4 + 39408 r 8 - 1 1849 r 9 

31104 r 3 

where A(0>(Gi, Mi, L 2 ,Qi, N 3 , L 4 , L 5 , M 3 ,G 6 )) = - [^3(4 ^3r 2 i/ + 8 ^3i/a- + 4r 2 y 2 - 16 y- 4r 3 ^3y- 24 y 2 + 12 x + 
24 r- 6r 2 - 12 xV - 12 r 3 x + 3 r 4 x 2 + 12 r 2 i + 20 \fly + 2 r 4 \/3y x + r 4 y 2 - 24)] / /[24r 2 



Case 11: 

p(x 2 e N r PE (x 1 )nr r 1 (x 1 ),x 1 g t.) = 

S11 f eam(x) + f 3w r«W + /•'" pW\ ^(^(G 1 ,M 1 ,L 2 ,Q 1 ,Q 2 ,L 5 ,M 3 ,G 6 )) ^ ix _ 

; A{T{y 3 )y 

(r - 1) (1080 r 16 + 1080 r 15 - 17820 r 14 - 540 r 13 + 65394 r 12 - 46926 r 11 + 105435 r 10 - 261765 r 9 + 229286 r 8 - 180586 r 7 + 
101638 r 6 + 40774 r 5 - 46112 r 4 + 24448 r 3 - 20224 r 2 + 10496 r - 6144) j j [l0368 r 3 (2 r 2 + l) 3 ] 

where A(&>(Gi, Mi, L 2 , Qi, Q 2 , L B , M 3 , G 6 )) = - [^3(6x + 3r 2 -4r 2 a; ^/3y-4y 2 - 6iV 2 + 2r 4 ^3yx + 4 ^/3y:r-2r 2 y 2 - 

4 r 3 ^y + r 4 y 2 - 12 r 3 x + 3 r 4 x 2 + 12 r 2 x - 6 + 4 \/3V 2 y + 2 v^y)] / [l2 r 2 ] . 
Case 12: 

,1/2 r r 3 (x)^ 



fl/2 rr 

P(x 2 g iVj s (Jfi) n r{(Xi),Xi eT s )= / 



A(T(y 3 )y dydx = 



(49 r 6 - 204 r 5 + 476 r 4 - 768 r 3 - 8 r 2 + 768 r - 288) (—6 + 5i 



7776 r 2 

where A{9>{Gy, G 2 , Qi, JV 3 , £ 4 , £5, M 3 , G 6 )) = - [\/3(-12a; + 12 r - 3r 2 + 12 r x - 20 V3r y - 12 xV + 4 ^3r 2 y - 12 r 3 x + 

3 r V + 28 %/3y + 12 x 2 + 12 r 2 :r - 12 - 20 y 2 + 4 r 2 y 2 - 4 r 3 V^y + r 4 y 2 + 2 r 4 \/3y a;)] / [24 r 2 ] . 
Case 13: 

(4r 7 + 8r 6 - 37 r 5 - 58r 4 -84r 3 + 168 r 2 +336r - 352) (-2 + r) (r 2 + 2r - 4) 2 

384 (r + 2) 2 r 2 

where A(0>(Li, L 2 , Qi, N 3 , L 4 , L 5 , L e )) = - [^(-4r 3 \/3y- 8 \/3r y + 12 z + 24r- 8y/3y x- 12 r 2 + 24r x -24- 12 x 2 r 2 + 

4 V3r 2 y - 32 y 2 - 12 r 3 x + 3 r 4 ^ 2 + 20 \/3y - 24 x 2 + 12 r 2 x + 2 r 4 ^/3y a + r 4 y 2 + 4 r 2 y 2 ) j |^24 r 2 J . 
Case 14: 

\J*11 Jr 12 (x) Js 14 Jr 10 (x) J A U \y3)) 

- [(135 r 11 + 675 r 10 - 1350 r' J - 9450 r 8 + 702 r + 39150 r 6 + 24272 r 5 - 47432 r 4 - 135040 r 3 + 57088 r 2 + 204800 r- 



134144) (r - 1)] / [l0368 (r + 2) 

where A(,^(Li, L 2j Qi, Q 2 , -^5, Le)) = - [V3(-4 r 3 V3y + 4 V3r y + r 4 y 2 + 6x-4: V3y x + 2 r 4 V3y x + 12rx-4r 2 x V$y 
6 x 2 r 2 + 4 V3r 2 y - 12 r 3 x + 3 r 4 a; 2 + 2 ^3y - 12 x 2 + 12 r 2 x - 6 - 8 y 2 - 2 r 2 y 2 )] / [l2 



Adding up the P(X 2 G Np E (Xi) D ri(Xi),Xi £ T,) values in the 14 possible cases above, and multiplying by 6 we 
get for r G [1,4/3), 

(r - 1) (5 r 5 - 148 r 4 + 245 r 3 - 178 r 2 - 232 r + 128) 
MandW = 54 r 2 (r + 2) (r + 1) ' 

The ^ an d(r) values for the other intervals can be calculated similarly. For r = oo, ^i an d(r) = 1 follows trivially. 



Derivation of iy an d(r) in Theorem 14.21 

By symmetry, P({X 2 ,X 3 } C JVJ a (Jfi) PI ri(Xi)) = 6P({X 2 ,X 3 } C AT^(Xi) n rftJfi), Ii£T s ). 

For r £ [6/5, yo — l), there are 14 cases to consider for calculation of f a nd(f) in the AND-underlying version: Case 1: 

P^XajcAW^nr^), x ie T.)= | r fj miX) + f° fj^ Am %T{y?)Y Gfi))2 dydx = 

(r 2 + l) 2 (r + l) 2 (r-l) 2 
384 r 10 

where A(0>(Gi,Ni,N a ,Ge)) = ^ {V^y + 3xf r 2 /36 - ^^j^ . 
Case 2: 

P({* 2 ,x 3} c ^(Xx) n iW), * e r.) = f f f W + f f ^ A( ^ (Gl ^' P2 'f' G6))2 ^ = 



s 5 Jr 5 (i) J 8 6 JO 



A(T(y 3 )) 3 



(5 + 38 r + 137 r 2 + 320 r 3 + 552 r 4 + 736 r 5 + 792 r 6 + 640 r 7 + 407 r 8 + 178 r 9 + 35 r 10 ) (-1 + rf 

960 r 10 (r + l) 5 

where A(£*(0?i, iVi, P 2 , M 3 , G 6 )) = - ^ ■ 

Case 3: 

p({x 2 ,x :i } c JV5 B (xo nrftxo, Xi e t.) = 

s 5 Jr 7 (x) Js 9 Jo Js 12 Jo J A{T(y S )) 3 

17496 r la - 122472 r 18 + 139968 r 17 + 524880 r 16 - 553095 r 15 - 595971 r 14 + 368826 r 13 - 724758 r 12 - 543876 r 11 + 
1416996 r 10 + 1646470 r 9 + 92870 r 8 + 523048 r 7 - 768368 r 6 - 1729902 r 5 - 1434990 r 4 + 122185 r a + 941941 r 2 + 

573440 r + 114688] / [2099520 (r + l) 5 r 10 ] 

where A(3 s (Gi,G 2 ,Qi, P 2 , M 3 ,G 6 )) = - [^5(4 v / 3r 2 y - 8r 3 v 7 ^ + 4r 2 y 2 + 4r 4 y 2 + 4j/ 2 + 8rV3j/a; + 6 - 12 xV - 12 x - 

12 r - 24 r 3 a; + 12 r 4 a; 2 + 9 r 2 + 12 r x - 4 v^r 2/ + 12 x 2 + 4 + 12 r 2 a;)] / [24 r 2 ] . 
Case 4: 

p({x 2 ,x 3 } c iv5 B (Xi)nrl(Xi), Xi e t.) = 

rraW^ r.io + /••«, Aj&jGi, M u L 2 , Q U P 2 , M 3 , G 6 )) 2 _ 

-'sg Jrg(x) J s 5 J r s (x) Ja 10 J r 3 (x) 

- [32768 - 409264128 r 7 + 1455989508 r 12 + 680709729 r 8 - 4423680 r 3 + 155509 r 2 + 22889801 r 4 + 202936917 r° + 
6011901 r 20 + 1060982949 r 16 - 614739456 r 17 + 240330993 r 18 - 56097792 r 19 - 77783040 r 5 - 999857664 r 9 + 
1299257316 r 10 - 1461851136 r 11 - 1407624192 r 13 + 1414729905 r 14 - 1352392704 r 15 ] / [2099520 (r 2 + l) 5 r 10 ] 

where A(&>(Gi, Mi, L 2 , Qi, P 2 , M 3 , G 6 )) = -[\/3(-6i- 2 r 2 - 3+6z- 12r 3 x + 6r 4 :r 2 - 4r 3 v / % + 4 V3yi + 4r 4 v / 3y2: + 

2 r V + 3 r 2 + 2 v 7 ^ 2 ?/ - 2 v^J/ + 2 7~V + 6 r- 2 z)] / [l2 r 2 ] . 
Case 5: 

P((X X\rN r QMnmi X cr\ ( H r (X) i r r 8W V(^(Gi^i,A,P 2 ,M 3 ,G 6 )) 2 

[(35361 r 16 - 229392 r 15 + 602820 r 14 - 858384 r 13 + 778848 r 12 - 460368 r 11 + 277740 r 10 - 258768 r 9 + 160594 r 8 - 62256 r- 
5892 r 6 - 17712 r 5 + 19224 r 4 + 11664 r 3 + 5076 r 2 + 1296 r + 405) (-12 r + 7r 2 + 3) 2 ] / [699840 r 10 (r 2 + if 

l. AfriBfsl 1/ n n i, ^^ ^(-4r 3 ^j-12r 3 1 +3 r 2 +6 r 4 ^ l/I 4-9r 4 I 2 4-3 r 4 ! , 2 f l , 2 +2V3 SI +3x 2 ) 

where A{&>{Gi, Mi, Pi, P 2 , M 3 , G 6 )) = 1272 i - 



Case 6: 



P({x 2 ,x 3 } c N£ E (Xi) n rl(Xi), Xi g t.) = 

/ rs 3 ri am (x) ^ ,. T + «,(.)\ A(^(Gi, Mi, Pi, P 2 , M 3 , G 6 )) 2 rf ^ _ 

W s 2 ■/ 1-5(1) 

- [3645 - 17496 r + 5003898912 r 28 + 31646646384 r 26 + 1100 9 8 9 44 r 30 - 1090803456 r 29 - 1463 75 1360 r 27 + 66339 r 2 - 
99072645696 r 23 + 79269457632 r 24 + 66 73 158 r 8 - 4870743552 r 13 - 168073488 r 9 + 535086 r 4 - 262440 r 3 - 1737936 r 5 - 

18592416 r 7 - 107383563504 r 21 - 41219053272 r 17 + 58981892347 r 18 - 78265758888 r 19 + 95887286866 r 20 + 
109053166552 r 22 + 5500 548 r 6 + 466565130 r 10 - 1070573040 r 11 + 2380992104 r 12 + 9191633420 r 14 - 16312513248 r 15 + 

26801184917 r 16 - 54759787776 r 25 J / [l399680 (r 2 + l) 5 (2r 2 + l)V°] 

where ^(^(Gi, Mi, Pi, P 2 , M 3 , G 6 )J = 1 12733 -• 

Case 7: 

p({x 2 ,x 3 } c JV5 B (xo niT(Xi), Xi e t.) = 

/ fS8 r9 (x) fS10 r°W\ A{&(G 1 ,M 1 ,L 2 ,Qi,P 2 ,M 3 ,G 6 ))\ , 

wow dydx = 



Jrg(x) J s 8 J r 2 (x) 



[4 (162576 r 22 - 1083456 r 21 + 3368016 r 20 - 69 6 9 8 88 r 19 + 11578088 r 18 - 15664080 r 17 + 18796852 r 16 - 19984824 r 15 + 
19534445 r 14 - 18170472 r 13 + 15507752 r 12 - 13150464 r 11 +9987958 r 10 - 7448736 r 9 +5016464 r 8 - 2991768 r 7 + 1857485 r 6 - 
749160 r 5 + 481804 r 4 - 96720 r 3 + 76160 r 2 - 4032 r + 4320) (2r - l) 2 (r - l) 2 ]/ [32805 (r 2 + l) 5 r 6 (2r 2 + l) 5 ] 

where A(0>(Gi, Mi, La, Qi, P2, M 3 , Ge)) = -[V3(-6a:V - 3 + 6z- 12r 3 :r + 6r 4 x 2 - 4r 3 ^y + 4\/3yx- + 4r 4 ^3ya; + 

2 r V + 3 r 2 + 2 %/3r 2 ?/ - 2 v/3j/ + 2 r- 2 y 2 + 6 r 2 x)] j [l2 r 2 ] . 
Case 8: 

p({x 2 ,x 3 } c jv£ B (Xi)nrT(Xi), Xi e t.) = 

/ rrsW^ ri/2 r a W\ A(^(Gi, G 2 , Qi, jV 3 , M c , M 3 , G 6 )) 2 d ^ = 

- [-458752 + 811008 r 2 + 329205504 r S - 582626304 r 13 - 489563136 r 9 - 65536 r 4 - 168708096 r - 57883680 r 17 + 
18009258 r 18 - 3623400 r 19 + 352563 r 20 + 4150 2 720 r 6 + 659111904 r 10 - 761846400 r 11 + 725173376 r 12 + 409477188 r 14 - 

254829600 r 15 + 135968852 r 16 ] / [8398080 1 



where A( <^(Gi, Ga, Qi, JV 3 , Mo, M 3 , G 6 )) = -[>/3(-12arV - 12 x - 12 r - 12r 3 a; + 3r 4 a; 2 + 4^3r 2 y + 5r 2 + 12 r a; + 

12 x 2 + 2 r 4 ^?/ x + 4 rV - 4 r 3 \/3y + 6 + 4 y 2 + r 4 y 2 + 4 \/3y + 12 r 2 x- - 4 \/3r y) j j [24 r 2 ] . 
Case 9: 



.6 T.) = 



p({x 2 ,x 3 } c iv; B (Xi)nrr(Xi), Xi . 

/ ~„ «(-) p «(-)\ A(^(Gi, Afi, L 2 Qi iV 3 M C , M 3 , G 6 )) 2 = r ia _ 49392 15 

170226 r 14 - 392112 r 13 + 680784 r 12 - 1040256 r 11 + 1385628 r 10 - 1337760 r 9 + 816224 r 8 - 253824 r 7 + 
469088 r 6 - 1029888 r 5 + 820992 r 4 - 488448 r 3 + 190976 r 2 + 49152 r + 8192) (- 12 r + 7 r 2 + 4) 2 ] / [8398080 r 10 ] 

where A(&>(Gi, Mi, L 2 , Qi, N 3 , M c , M 3 , G 6 )) = - [\/3(-12arV - 6 + 12 x- 12r 3 x + 3r 4 x 2 + 2r 2 + 2r 4 V3yx + r 4 y 2 
8 V3y x + 4 r 2 y 2 - 4 + 4 \/3r 2 y + 12 r 2 x - 4 r 3 \/3y)] / [24 r 2 ] . 



Case 10: 

p{{x 2 ,x 3 } c ^(JfOnr^ij), Xi e t.) = 

/ r.i4 pW + /-M pW + r 1 / 2 /ri2(-)\ A(^(Gi, Mi, L 2 , Qi, JV 3 , £ 4 , £5, M 3 , G 6 )) 2 rf ^ = 

W*10 J r 2 (x) J S14 Jr 2 (x) J s 13 J r 3 (x) 

[4423680 - 4627454976 r 6 + 511684992 r 11 + 2163142656 r 7 - 660127744 r 2 - 31555584 r + 3534520320 r 3 + 7647989760 
7785504 r 15 - 1313880 r 16 + 19683 r 18 - 7240624128 r 4 - 1511047552 r 8 + 1204122240 r 9 - 796453824 r 10 - 282583320 ? 

107804736 r 13 - 30362052 r 14 ] / [l6796160 r 6 ] 



8 1 
r + 

12 . 
r + 



here A(^ , (Gi,M 1 ,L2,Qi,N 3 ,L 4 ,L 5 ,M 3 ,G( i )) = - U/z{-l&V3ry + 20V3y - 24 y 2 - 12x 2 r 2 + 12 x + 24r - 12r 3 a 

rV - 6 r 2 - 24 + 4 VSr 2 y + 8 V3y x-4 r 3 V3y + 4 r 2 y 2 + r A y 2 + 2 r 4 ^y x + 12 r 2 x)] / [24 r 2 ] . 
!ase 11: 

p({x 2 ,x 3 } c iv5 B (Xi)nrl(jri), Xi e t.) = 

/ rW«0 rio rna(») r. M pW\ A(g»(Gi, Afi, L 2 , Qi, Q 2 , L 5 , M 3 , 

Ws 7 J r (x) Jen J r 9 (x) J s 10 J r 10 (x) 



(r - 1) (-1474560 + 8847360 r + 111456 r 26 + 111456 r 27 - 27738112 r 2 + 23311152 r 23 - 167184 r 24 - 808889416 r 8 - 
2228253688 r 13 + 366739256 r 9 - 207619072 r 4 + 98557952 r 3 + 397199360 r 5 + 802401664 r 7 - 34733448 r 21 - 624736557 r 17 + 
400615470 r 18 - 134938386 r 19 + 39014136 r 20 - 180 2 6 64 r 22 - 6400 5 8 432 r 6 + 407655352 r 10 - 1227078728 r 14 + 
1996721576 r 12 + 2033409092 r 14 - 1681870468 r 15 + 1064030499 r 16 - 2842128 r 25 )] / [l866240 (2r 2 + l) 5 r 6 ] 

where A(^(Gi, Mi, L 2 , Qi, Q 2 , £5, M 3 , G 6 )) = - \S(4 V3r 2 y + 4 V$y x - 2r 2 y 2 - 4r 3 V3y - 4y 2 - 1 V3r 2 yx - 6x 2 r 2 + 

6 x - 12 r z x + 3 r 4 x 2 + 3 r 2 + 2 r 4 \/3jf a: + r 4 y 2 + 2 \/3y + 12 r 2 x - 6)] / [l2 r 2 ] . 
Case 12: 

[(2322432 - 7554816 r + 9510912 r 2 + 1046068 r 8 - 558720 r 9 + 2444224 r 4 - 5799360 r 3 - 2134656 r 5 - 1608672 r 7 + 

2169696 r 6 + 216300 r 10 - 55440 r 11 + 7095 r 12 ) (-6 + 5 r) 2 ] j [4199040 r 4 ] 

where A(^»(Gi, G 2 , Qi, JV 3 , L 4 , L 6 , M 3 , G 6 )) = - [V3~(-12a;V - 12 a; + 12 r - 12r 3 x + 3r 4 x 2 - 3r 2 + 12 r x + 28 ^3y + 

12 x 2 - 20 1/ 2 + 12r 2 x + rV +4r 2 y 2 - 4 r 3 V3y + 2 r 4 V3y x + 4 V3r 2 y - 20 y/Sry - 12) j j [24 r 2 ] . 
Case 13: 

y si4 JruW ^(T(^3)) 3 

- [(9r 14 + 36 r 13 - 132 r 12 - 576 r 11 + 164 r 10 + 2512 r 9 + 4976 r 8 - 1536 r 7 - 13888 r 6 - 17536 r 5 - 3072 r 4 + 79360 r 3 + 

9216 r 2 - 120832 r + 61440) (-2 + r) (r 2 + 2 r - 4) 2 ] j [768O (r + 2) 3 r 4 ] 

where A(9>(Lx, L 2 , Qi, iV 3 , L 4 , L B , ie)) = - [^3(r 4 y 2 -8^3ry-8^3j/x + 4r 2 y 2 -4r 3 ^3y-32y 2 + 2r 4 ^3i/j--12x 2 r 2 + 

12 x + 24 r - 12 r 3 x + 3 r 4 x 2 - 12 r 2 + 4 ^r 2 y + 24 r x - 24 a; 2 - 24 + 20 V%y + 12 r 2 x)^ j [24 r 
Case 14: 



p{{x 2 ,x 3 } c JV5 B (Xi)nri(Jsri), Xi e r.) = 

/ ,. 14 r/.„(«) ,1/2 ^. m (-)\ A(^(L 1 ,L 2 ,Qi,Q a ,L 5 ,L 6 )) 2 J _ 

[(r - 1) (3483 r 18 + 24381 r 17 - 34830 r 16 - 529416 r 15 - 265680 r 14 + 4274208 r 13 + 4999320 r 12 - 15227352 r 11 - 
25751336 r 10 + 19466488 r 9 + 62834064 r 8 + 17452256 r 7 - 53339200 r 6 - 117114624 r 5 - 51206656 r 4 + 270430208 r 3 + 

58073088 r 2 - 296222720 r + 122159104) ] j [l866240 (r + 2) 3 r 4 ] 



where A{@>(L X , L 2 , Qi, Q 2 , L 5 , L 6 )) = -[y/3(-4y/3yx -2r 2 y 2 + 4^3ry - 4r 3 ^3y - 8y 2 ~Ay/?>r 2 yx - 6x 2 r 2 + 6x- 
12 r 3 x + 3 r 4 x 2 + 4 V3r 2 y + 12 r x - 12 x 2 + 2 r 4 v^J/ x + r V + 2 v 7 ^ + 12 r 2 x - 6)] / [l2 r 2 ] . 

Adding up the P({X 2 ,X3} C JVp B (-Xi) l~l ^i(Xi), Xi G T s ) values in the 14 possible cases above, and multiplying by 
6 we get for r£ [6/5, y/E — l) , 

iw(r) = - [219936 r - 3041936 r 2 - 30889822 r 8 + 18084672 r 13 + 27137438 r 9 + 2364868 r 4 + 2305864 r 3 - 4168820 r 5 - 
2832544 r +486 r 21 - 118850 r 17 - 45155 r 18 - 269 r 19 + 3402 r 20 + 11101160 r 6 + 24604048 r 10 - 43009544 r 11 +8770788 r 12 - 

13736295 r 14 + 2751855 r 15 + 443518 r 16 + 49152] j [ll6640 r 6 (r + 2) 2 (2 r 2 + l) (r + l) 3 ] . 

The Va.ndir) values for the other intervals can be calculated similarly. 

Appendix 3.2: Derivation of n OT (r) and ^ or (r) for Uniform Data 

Derivation of fi or {r) in Theorem 14.21 

First we find /tior(f) for r G [l, oo). Observe that, by symmetry, 

M or(r) = P(X 2 G JV^ B (Xi)urj(Xi)) = 6P(X 2 G A?5(Xi)url(Xi),Xi G T s ). 

For r G [1,4/3), there are 17 cases to consider for calculation of v OT (r) in the OR-underlying version. Each Case j 
correspond to Ri for i = 1, 2, . . . , 17 in Figure [T6l Case 1: 

p( y a N r ( Y \ \ \ ~r r ( Y } x a t \ ( H [ lam{X) + f 1 f eamix) \ A{&{A,My,M c ,M 3 )) 



± r 2 _4r/9 + l/3 



where A(£P(A, Mi, Mc , M3)) = v^/12. 
Case 2: 

P(X 2 G N£ E (Xy) U rl'(Xi),Xi G T s ) = 



ri(t) + p /"■»<*> + /•«* r* w + T 5 r 5( *°] AQg(A, m x ,l 2 ,l 3 , m c , m 3 )) _ 

JO Jsi JO " S3 JO 

(r - 1) (1817 r 7 - 7807 r 6 + 14157 r 5 - 14067 r 4 + 7893 r 3 - 2475 r 2 +405 r - 27) 

864 r 6 

V^3(-4 r + 12n + 5 r 2 +3y 2 +6 \/3!/— 6 \/3jl 21+9-18 a; + 9 a: 2 ) 



where My, L 2 , L 3 , M c , M 3 )) = 

Case 3: 



P(Y a N r IY \ \\ T r ( Y \ Y a T \ ( H J, H f A{PJ>{A,G 2 ,G 3 , M 2 , M c , M 3 )) 

P(X 2 G NMXy) U F l( Xi),Xi G T.) = ^ ^ +J s5 ^ j j^-^ dydx = 

(I3r 4 - 4r 3 + 4r - 1 - 2r 2 ) (r - l) 4 



96 r 6 

where A(&>(A, G 2 , G 3 , M 2 , M a , M 3 )) = * ^ -■ 

Case 4: 

f f 2 p"^ 3 , f 3 f 5 ^ ^(AMi,l 2 ,l 3 ,L4,l 5 , Af 3 )) , 

P(X 2 G A^(Xi) U r l( Xx), * G T.) = +y^ j wm - 2 dydx = 

(9 - 72 r + 192 r 2 - 192 r 3 + 76 r 4 ) (4 r - 3 + V^f (4 r - 3 - \/3) 2 

10368 r 6 

u ^/^i/^ AT r T T r 1M\ 73(4 % /3r ! /+9r 2 - 2 4r+l 2 r 3: +15H 2 -6^3!/-6^3 ! ; 3: +18-18 a :+9 a: 2 ) 

where A{&>{A, My , L 2 , L 3 , L 4 , L 5 , M 3 )) = — i . '-. 



12 r 



Case 5: 

P (Y <= ^ miim \ x a t \ ( r r (x \ r r (x) \ m&(ag 2 ,g 3 ,m 2 ,m c ,p 2 ,n 2 )) 

p(x 2 e n pe (Xi) u r^),* g r.) = ^ +y s8 ^ j dyrf* = 

(-1 + 2r + 6r 2 - 6r 3 + 22r 5 + 17 r 6 ) (r - l) 3 
96 r 6 (r + l) 3 

where A(0>(A, G 2 , G 3 , M 2 , M c , P 2 , N 2 )) = [^3(-2y 2 - 4^3y + 4^3yx - 6 + 12 x -6x 2 + 7r 2 - 4r 3 ^3y - 12r 3 x + 

8 r 4 V3y x + 12 r 4 x 2 + 4 r 4 y 2 ) j / [24 r 2 ] . 
Case 6: 

P(X 2 6 A^PG) U rl(Jfi),Xi G T.) = 

/ p raW^ « u ^ ri/a ,r 8 (s)\ ^(A, iVi,Qi, G 3 ,M 2 ,M C ,P 2 , N 2 )) _ 

; ^(T(^ 3 )) 2 

81 r 9 - 189 r 8 + 561 r 7 - 45 r 6 - 1894 r 5 - 18 r 4 + 1912 r 3 + 224 r 2 - 384 r - 128 

1296 (r + l) 3 r 4 

where A(^(A, iVi , Qi, G 3 , M 2 , Mc, P>, AT 2 )) = [^3(4r y 2 - 4 \/3y + 12 x + 13 r - 12 + 18 r 3 a; 2 + 12 r x - 12 r a; 2 - 8 \/3r 2 y + 

4 \/3r y-24r 2 x + 12 V3r 3 y x + 6 r 3 y 2 )] / [24 r] . 
Case 7: 

p(x 2 e iv; B (jfi) u r[(Xi),Xi g r.) = 

/ p + po rmCO + /-.la A(g»(A, N 1 ,Q 1 ,L 3 , Mg, P ^ N ^)) dydx = 

- [l28 - 1536 r - 302592 r 7 + 11753 r 12 + 346171 r 8 - 28416 r 3 + 8384 r 2 + 69760 r 4 + 220201 r 6 - 135936 r 5 - 305664 r 9 + 

186683 r 10 - 69120 r 11 ] / [l944 (r 2 + l) 3 r 6 ] 

where A(&(A, Ni, Qi, L 3 , Mc, P 2 , iV 2 )) = [^3(-4 y + 2 \/3> 2 y - 12 a; - 12 r + 8r 2 + 12 r a; - 6a; 2 r 2 + 2 r 2 y 2 - 4 ^3y a- + 

3r 4 y 2 -4r 3 ^3y- 12r 3 x + 9r 4 x 2 + 4 ^3y + 6 r 4 ^3y x + 6 a; 2 + 6 r 2 x + 6 + 2 y 2 )] j [l2 r 2 ] . 
Case 8: 

P(* 2 G W ) L) T\{X\),X\ € T s ) = ( f" f W + f ^(^(A iV, A y 3 Me, P 2 , iV 2 )) = 

Visa y S8 JrsM y -4(r(^3)) 2 

(895 r 8 - 2472 r + 3363 r 6 - 2880 r 5 + 2220 r 4 - 1296 r 3 + 675 r 2 - 216 r + 27) (- 12 r + 7 r 2 + 3) 2 

7776 (r 2 + l) 3 r 6 

where A(J*(A, JVi, Pi, L 2 , L 3 , M C , P 2 , N 3 )) = [^3(4r 4 y 2 + 8r 4 ^3yx + 12r 4 a; 2 -4r 3 ^3y-12r 3 a;-4^3ry-12r + 12ra; + 

8 r 2 + 3 y 2 + 6 ^3y - 6 ^3y x + 9 - 18 x + 9 a; 2 ) j / [l2 r 2 ] . 
Case 9: 

P(X 2 G iVpBpG) U ri(Xi),JCi G T s ) = 

/ r. a /•<.»{«) + r.T ri-mW + r.« A(^(A, JVi, P 1; L 2 ,L 3 ,L 4 , L 5 , P 2 , N 2 )) = 

\ J s 2 J r 5 (x) J s 3 J r 2 (x) •* s 7 J r 2 (x) 

/ ^(T(^3)) 2 

- [355328 r 18 - 2204160 r 17 + 6591792 r 16 - 13254912 r 15 + 20639832 r 14 - 26417664 r 13 + 28578916 r 12 - 26760576 r n + 
21960774 r 10 - 15877152 r 9 + 10180620 r 8 - 5753232 r 7 + 2856483 r 6 - 1222128 r 5 + 438777 r 4 - 128304 r 3 + 28107 r 2 - 

3888 r + 243] / [7776 (r 2 + if (2 r 2 + l) 3 r 6 ] 

where A(£?(A, N ly P lt L 2 , L 3 , L 4 , L 5 , P 2 , N 2 )) = k/3(l8 + 4^3ry - 18a; - 24r + 12 r 2 + 12rx - 6^3y + 8r 4 ^ya; - 
12 r 3 a; + 12 r 4 a; 2 +9x 2 + 15y 2 +4 r 4 y 2 - 4 r 3 ^3y - 6 \/3y sc)] j [l2 r 2 j . 



Case 10: 

p(x 2 e n pe (x.) u w),* e r.) = ^ +y^ y^ j dydx = 

! r 12 - 864 r 11 + I486 r 10 - 1896 r 9 + 2056 r 8 - 1608 r 7 + 1189 r e - 654 r 5 + 317 r 4 - 132 r 3 + 44 r 2 - 12 r + 2) (2 r - l) 2 

(r-l) 2 ]/[243 (r 2 + l) 3 (2r 2 + l)V] 

where A(^>(A, JVi, Qi, L 3 , U, L 5 , P 2 , JV 2 )) = [\/3(4 %/3r y + 2 v^r 2 ?/ - 8 \/3j/ - 12 x - 24 r + 12 r 2 + 12 r x - 6 xV + 15 - 
12 r 3 x + 9 r 4 x 2 + 6 x 2 + 6 r 2 x + 6 r 4 ^3y x + 2 r 2 j/ 2 - 4 \/3j/ x + 3 r 4 y 2 - 4 r 3 V3y + 14y 2 )jy / [l2 r 2 ] . 
Case 11: 

p ( x 2 e amxou gT . ) = ( r 3 r w + f 2 r w> ) .y,,*,,*,)) 

_ 1536 - 6528 r 2 + 133834 r 8 - 48240 r 9 + 95616 r 4 - 20736 r 3 - 158976 r 5 - 200064 r 7 + 196680 r 6 + 7107 r 10 

15552 r 4 

where A{&{A, iVi, Qi, G 3 , M 2 , N 3 , JV 2 )) = [^3(4 ry 2 + 12x + 9r-12 + 9r 3 x 2 + 12 rx - 12 rx 2 - 4^3r 2 y + 4^3ry + 

6 \/3r 3 y x + 3 r 3 y 2 - 12 r 2 x - 4 y/3y)^ j [24 rj . 
Case 12: 

(147 r 8 - 504 r 7 + 530 r b - 336 r 5 + 876 r 4 - 1056 r 3 + 896 r 2 - 384 r + 64) (- 12 r + 7 r 2 + 4) 2 

15552 r 6 

where iVi, Qi, L 3 , iV 3 , N 2 )) = [^3(4y 2 - 8^3yx - 24x - 24r + 8^% + 12 r 2 + 4^3r 2 y + 6r 4 ^3yx + 24rx - 

4r 3 V3y + 3r 4 y 2 - 8 V3r y - 12 x 2 r 2 - 12r 3 x + 9r 4 x 2 + 12x 2 + 12r 2 x + 4r 2 y 2 + 12)] j [24 r 2 ] . 
Case 13: 

p(x 2 e ^(xOur^Xx),^! e r.) = 

811 r ioW + r 13 r i2W + r 1/2 r^*A a(^(a,jvi,Qi,z,3,jv3,jv 2 )) = 

; A(r(3> 3 )) 2 

[l024 - 12288 r + 295680 r 7 + 1053 r 12 - 197140 r 8 + 626688 r 3 - 100864 r 2 - 1294848 r 4 - 686528 r 6 + 1282560 r 5 + 

114336 r 9 - 30930 r 1() ] / [31104 r 4 ] 

where A(0>(A, Ni, Qi, L 3 , N 3 , N 2 )) = [^3(4y 2 - 8V3yx - 24x - 24 r + 8-^/3?/ + 12 r 2 + 4\/3r 2 j/ + 6r 4 ^?/x + 24 rx - 

4r 3 V3y + 3r 4 y 2 - 8 \/3V y - 12 a;V - 12r 3 x + 9r 4 x 2 + 12x 2 + 12r 2 x + 4r 2 ?/ 2 + 12) j / [24 r 2 j . 
Case 14: 

p(x 2 e ivj B (Jsri)ur;(Xi),Xi e r.) = 

ai1 + /••«> /••■«(«) + 4(^(A,jV 1 ,Q 1 ,L 3 ,L 4 ,Q 2 ,iV 2 )) rf ^_ 

; ^(r(^)) 2 

- [(r - 1) (1512 r 17 + 1512 r 16 - 16740 r 15 + 540 r 14 + 84078 r 13 - 83538 r 12 - 164835 r 11 + 401085 r 10 - 487872 r 9 + 
535728 r 8 - 463124 r + 335596 r 6 - 197440 r 5 + 64640 r 4 - 7936 r 3 - 1792 r 2 + 5632 r - 512)]/ [5184 (2r 2 + l) 3 ?" 4 ] 

where v4(^(A, JVi, Qi, L 3 , L 4 , Q 2 , iV 2 )) = [\/3(-6 x- 12 r + 6r 2 +6 r x + 2 ^r 2 y-r 2 ?/ 2 -2 V3yx+r 4 y 2 + 5 y 2 -2 r 2 x V3y + 

2 r 4 V3y x + 2 V3r y-2 r 3 V3y - 3 x 2 r 2 - 6 r 3 x + 3 r 4 x 2 - 2 \/3y + 3 x 2 + 6 r 2 x + 6)] / [6 
Case 15: 

.1/2 frsW a(^(A,JVi,Qi,G 3 ,M 2 ,JV 3 ,JV 2 )) 



P(x 2 g Arj s (Xi) u r;(jfi),Xi eT B )= / 

-V Si 3 «/ T* 1 



A(T(y s )r dydx - 

(147 r 5 - 612 r 4 + 980 r 3 - 768 r 2 + 744 r - 288) (-6 + 5 r) 2 
7776 r 



where A(9>{A, N u Qx, L 3 , L 4 , Q 2 , N 2 )) = [^3(4 ry 2 + 12x + 9r-12 + 9r 3 x 2 + 12rx - 12rx 2 - 4 ^flr 2 y + 4 ^3r j/ 
6 \/3r 3 ?/ s + 3 r 3 y 2 - 12 r 2 x - 4 v/3y)] / [24 j 



Case 16: 

P(x 2 g n pb (Xx) u r^),* e r.) = f 2 A( ^^;^ JV '- jVa)) d yd 3; 



; S14 Jru.v A{T{y z w 

(13 r 8 + 52r 7 + 10 r 6 - 184 r 5 + 60 r 4 + 624 r 3 - 48 r 2 - 832 r + 448) (-2 + r) (r 2 + 2r - 4) 2 

384 (r + 2) 3 r 2 

where A{8?{A, Nx, Qx, L 3 , N 3 , JV 2 )) = [^3(4y 2 - 8^3yx - 24x - 24r + 8 \/% + 12 r 2 + 4^3r 2 i/ + 6r 4 ^3yx + 24rx - 

4r 3 ^3y + 3r 4 y 2 - 8 \/3> j/ - 12 x 2 r 2 - 12 r 3 x + 9 r 4 x 2 + 12 a; 2 + 12r 2 x + 4r 2 y 2 + 12)]/[24r 2 ] . 
Case 17: 

P(X 2 6JV PS (Xi)uri(Xi),Xi 6T S )= / / +/ / ZT^rv dydx = 

\J S11 Jr 12 (x) Js 14 Jr 10 (x) ) y-yS)) 

[(189 r 12 + 1323 r 11 + 1026 r 1() - 10692 r 9 - 14364 r s + 51732 r 7 + 64664 r 6 - 183952 r 5 - 153504 r 4 + 398080 r 3 + 124928 r 2 - 

470528 r + 197632) (r - 1)] j [5184 r 2 (r + 2) 3 ] 

where A(^(A, Nx,Qx,L 3 ,N 3 , N 2 )) = [\/3(-6 x- 12 r + 6r 2 + 6r z + 2 \/3r 2 y - r 2 y 2 -2 V3y x + r 4 y 2 + 5 y 2 - 2 r 2 x v/% + 
2 r 4 \/3y x + 2 \/3r y-2 r 3 V3y - 3 xV - 6 r 3 x + 3 r 4 x 2 - 2 v^y + 3 x 2 + 6 r 2 x + 6)] / [6 r 2 j . 

Adding up the P(X 2 £ iVp B (Xi) U r'i(.X'i), .Xi G T s ) values in the 17 possible cases above, and multiplying by 6 we 
get for r € [1,4/3), 

_ 860r 4 - 195r 5 - 256 + 720r - 846r 3 - 108r 2 + 47r 6 
Vm V' ~ 108 r 2 (r + 2) (r + 1) ' 

The v OT {r) values for the other intervals can be calculated similarly. 

Derivation of u or (r) in Theorem 14.21 

By symmetry, P({X 2 ,X 3 } C N r PE (X 1 )ur r 1 (X 1 )) = 6P({X 2 ,X 3 } C N r PE (X 1 )ur r 1 (X 1 ), Xx € T s ). For r e [6/5,^5-1), 
there are 17 cases to consider for calculation of v OT (r) in the OR-underlying version (see also Figure [16]): Case 1: 

p({x 2 ,x 3} c n pb{ xx) u v U xx), xx e T .) = i r r (x) + r ^^^ ^ = 



s J rx(x) 



A(T(y 3 )r 



— r 2 r + 1/9 

81 27 1 



where A(0>(A, Mi, Mc, Ms)) = 1/12 \/3. 
Case 2: 

P({x 2 ,x 3 } c Af^fXOur;^), Xx <= t.) = 

pi pW p 3 pW pi pW p= pW\ A(,^M,Mi,L 2 ,L3,Afc,M 3 )) 2 , , 

/. I 7. ./, 7 3 ./ ' ./, 4 W J wmr ' /// ' / -' 

- [(r - 1) (119155 r 11 - 845345 r 10 + 2724777 r 9 - 5206743 r 8 + 6475257 r - 5454855 r 6 + 3155193 r 5 - 1249479 r 4 + 

332181 r 3 - 56619 r 2 + 5589 r - 243)] / [25920 r 10 ] 

where A(^»(A, Af a , L 2 , L 3 , M c , M 3 )j = — ^ ^ '-. 

Case 3: 

Pi{x 2 ,x 3} c n pe{ Xi) u FKxo, e t.) = (£ jT W + £ jT W ) ^%g^^ - 

(215 r 8 - 136 r 7 - 56 r 6 + 172 r 5 - 55 r 4 - 60 r z + 66 r 2 - 24 r + 3) (r - l) 4 

2880 r 1 " 

where A(^>(A, G 2 , G 3 , M 2 , M c , M 3 )) = K - -f-^ >-. 



Case 4: 



p({x 2 ,x 3 } c x£ B (Xi)url(Xi), Xi e t.) = 

S2 /-W») + p /"*<->\ A(^(A,M 1 ,L 2 ,L 3 ,L4,L5,M 3 )) 2 ^ d:r _ 

si ■/r 2 (a;) is 2 J ra(x) / -^(^O^)) 3 

[(37072 r 8 - 195072 r 7 + 453120 r 6 - 589248 r 5 + 460728 r 4 - 217728 r 3 + 60480 r 2 - 9072 r + 567) 

(4r-3 + V3y (4r-3- V^) 2 ] /[ 

u jr r r r r uH v^(4 \/3r y+9 r 2 -24 ^+12 r 1+15 y 2 -6 x/3i/-6 V3l/ ^+18-18 1+9 a: 2 ) 

where A{&>{A, Mi , L 2 , L 3 , La, is, M 3 )) = — => j^p -. 

Case 5: 

p({x 2 ,x :i } c ^(xour^Jd), x x g t.) = 

S6 r r 7W fS9 r r 7W \ A(^(A,G 2 ,G 3 ,M 2 ,Mc,P 2 ,X 2 )) 2 



1866240 r 10 



5 •''•5(2:) ^=6 "'0 



3 



dyda; = 



(3 - 12 r - 15 r 2 + 84 r 3 + 18 r 4 - 232 r 5 + 130 r 6 + 504 r 7 - 108 r 8 - 288 r 9 + 623 r 10 + 920 r 11 + 373 r 12 ) (r - l) 3 

2880 r 10 (r + 1) 5 

where A(&(A, G 2 , G 3 , M 2 , Mo, P 2 , X 2 )) = [%/3(-2y 2 -4^3*/ + 4V^yx - 6 + 12s - 6x 2 + 7 r 2 - 4r 3 V3y - 12r 3 x + 
8r 4 V3yx + 12r 4 x 2 + 4r 4 j/ 2 )] / [24r 2 l 
Case 6: 



G T 3 ) = 



p({x 2 ,x 3 } c JV5 B (Xi)urI(jfi), 

/ /••» /••« pW rV2 pW\ A(^(A,7V 1 ,Q 1 ,G 3 ,M2,M c ,P2,iV2)) 2 , , 

[L + L i + L I ) w dydx 

- [l9683 r 15 - 59049 r 14 + 83106 r 13 + 167670 r 12 - 211626 r 11 + 344466 r 10 - 142614 r' J - 2573586 r 8 - 128853 r 7 + 
3465675 r 6 + 1103824 r 5 - 1473304 r 4 - 730880 r 3 + 107776 r 2 + 158720 r + 31744] / [l049760 {r + l) 5 r 6 ] 

where A(J»(A, Ni, Qi, G 3 , M 2 , M c , P 2 , X 2 )) = [^3(4ry 2 + 12x + 13r + 12ra;-4^3y-12 + 4 V3r y - 8 ^3r 2 y + 18 xV - 

12 r x 2 + 6 r 3 y 2 - 24 r 2 s + 12 \/3r 3 y x)] j [24 r] . 
Case 7: 

P({X 2 ,X 3 } C iV5 B (Xi)un(Xi), Xi G T.) = 

/ rr,(-) + /-no p(-) + /••« Xi, Qi , £3, Af C , ft, X 2 )) 2 _ 

Wsg Jrg(x) J s 5 Jr s (x) J s 10 J r s (x) y ^(p(^)) 3 
- [6144 - 110592 r - 310846464 r + 2127553557 r 12 + 570050560 r 8 - 5031936 r 3 + 936960 r 2 + 19526656 r 4 + 147203072 r 6 + 
7627473 r 20 + 14190 720 42 r 16 - 762467328 r 17 + 288811029 r 18 - 68327424 r 19 - 59166720 r 5 - 923627520 r 9 + 
1340817105 r 10 - 1765251072 r 11 - 2350015488 r 13 + 2339575338 r 14 - 2016377856 r 15 j j [262440 (r 2 + l) 5 r 10 ] 

where A(^(A, Xi , Qi , L 3 , Mo, P 2 , X 2 )) = [\/3(-4 \/3r y + 2 \/3r 2 y - 61V - 12 x - 12 r - 12r 3 i + 9r 4 x 2 + 8r 2 + 12 r x + 

6 x 2 + 6 r 4 \/3y x + 2 r 2 t/ 2 - 4 \/3y x + 3 r 4 t/ 2 - 4 r 3 \/3?/ + 4 \/3?/ + 2 ?/ 2 + 6 r 2 x + 6)] / [l2 r 2 ] . 
Case 8: 

P({X 2 ,X 3 } C X; E (X 1 )un(X 1 ), Xi G T s ) = 

( r s , r 5 r sW \ a(^(a,xi,pi,l 2 ,l 3 ,m c ,p2,x2)) 2 
U 3 X 5W 4 Jr^J A(T(y 3 )r " ! " h ' 

[(426497 r 16 - 2443992 r 15 + 6726107 r 14 - 11753232 r 13 + 15220771 r 12 - 16367448 r 11 + 15754449 r 10 - 13773024 r 9 
10839672 r 8 - 7552440 r + 4592889 r 6 - 2374272 r 5 + 1018899 r 4 - 344088 r 3 + 81891 r 2 - 11664 r + 729) 

(-12 r + 7r 2 + 3) 2 ] / [699840 (r 2 + l) 5 r 10 ] 



r 9 + 



where A{S»{A, Ni, Pi, L 2 , L 3 , M c , P 2 , N 2 )) = [^(-4r 3 ^3y - 12 
12 r x + 3 y 2 + 6 V3y - 6 V3y x + 8 r 2 + 9 - 18 x + 9 x 2 )] j [l2 r 2 ] . 



r 3 x + 8r 4 ^yx + 12r 4 x 2 + 4r 4 y 2 - 4^3ry - 12 r + 



Case 9 

p({x 2 ,x 3 } c jv5 B (Xi)url(jfi), Xi g t.) = 

s 2 Jr 5 (x) J S3 Jr 2 (x) J s 7 J r 2 (x) 



33 + p p«(*> + p<*>\ A(g%A, JVx, A, L 2 ,L 3 ,L 4 ,L 5 ,P 2 ,N 2 )) 2 dydx _ 



- [l5309 - 367416 r + 60475010560 r 28 + 437704472832 r 26 + 1444872192 r 30 - 132 50 10 1248 r 29 - 18590 9 8 705 92 
4148739 r 2 - 2027754648576 r 23 + 1397612375040 r 24 + 20 429 1 77589 r 8 - 677278256112 r 13 - 49656902904 
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r + 



r 9 + 



)r 1T + 
)r i5 + 



48739 r - 2027754648576 r- + 1397612375040 r + 20429177589 r - 677278256112 r - 4965690291 
159963012 r 4 - 30005640 r 3 - 681714144 r 5 - 7515142416 r - 3097406755584 r 21 - 2609245249920 r 1 
3051035360256 r 18 - 3315184235136 r 19 + 3337272236928 r 20 + 263 194150 7968 r 22 + 24 3 5 97 1 806 r 6 - 
109069315047 r 10 - 218273842152 r 11 + 400534503738 r 12 + 1059615993384 r 14 - 1538314485120 r 15 

2076627064432 r 16 - 845838600192 r 25 ] / [l399680 (r 2 + l) S (2 r 2 + l) S r 10 

where A(9>(A, Ni, Pi, L 2 , L 3 , U, L 5 , P 2 , N 2 )) = [\/3(l8- 18 x- 24 r - 12 r 3 x + 12 r 4 x 2 + 12 r 2 + 12r a; + 4 v / 3ry-4r- 3 v / 3y 

4r 4 y 2 -6^3yx + 8r 4 V3yx + 9x 2 + 15 y 2 - 6 ^3y)j / [l2 r 2 ] . 
Case 10: 

p({x 2 ,x 3 } c N r PE {Xi) u rj(*i), Xi g t.) = 

/ p p W + p p W \ A(^(A,M,Qi,L, ^ L B ,ft,^))' 

[64 (12 - 144 r + 924 r 2 - 683328 r 23 + 1129 76 r 24 + 75 7211 r 8 - 10554918 r 13 - 1513230 r 9 + 16242 r 4 - 4320 r 3 - 51372 r 5 - 
344988 r 7 - 4867848 r 21 - 18583080 r 17 + 16493828 r 18 - 12883116 r 19 + 8668124 r 20 + 2 1 77536 r 22 + 14 1366 r 6 + 2774371 r 10 - 
4692510 r n +7331714r 12 +14002613r 14 -16948218r 15 +18708475r 16 ) (r - l) 2 (2r- l) 2 ]/ [32805 (r 2 + l) 5 (2r 2 + l) 5 r 8 ] 

where A(&>(A, Ni, Qi, L 3 , L 4 , L s , P 2 , N 2 )) = [\/3(2 V3r 2 y + 15 - 6xV - 12 x- 24r- 12 r 3 x + 9r 4 x 2 + 12 r 2 + 12 rx - 

8 V3y + 6 x 2 + 6 r 4 \/3y x + 14 y 2 - 4 \/3y x + 2 r a y 2 - 4 r s V3y + 3 r 4 y 2 + 6 r 2 x + 4 V^r y)] / [l2 
Case 11: 

P{{X 2 ,X 3 } C iVpB(Xi)ur^Xi), Xi G T s ) = 

/ /■sia rr 3 (x) t-1/2 rr 

\J 812 J r 6 (x) Js 13 Jr e {x) / 

- [-253952 + 1529856 r 2 + 601574256 r 8 - 385780320 r 13 - 776518272 r 9 + 7803648 r 4 - 70917120 r 5 - 396524160 r 
209710080 r 6 +869661288 r 10 -845940960 r 11 +668092108 r 12 + 147067614 r 14 -32610600 r 15 +3173067 r 16 ] / [8398080 r 6 j 

where A{&>(A,Ni,Qi,G 3 ,M 2 ,N 3 ,N 2 )) = [^3(4ry 2 + 12x + 4^ry + 9r - 4^3y + 12 r x - 12 + 9 x 2 r 3 + 6 V3r 3 y x - 

12 rx 2 -4V3r 2 y- 12 r 2 x + 3 r 3 y 2 )] /[24rj. 
Case 12: 

p({x 2 ,x 3 } c jv£ B (Xi) url(Xi), Xi g t.) = 

/ p « W+ p pw\ W y,y„^ 

[(64827 r 16 - 444528 r 15 + 1223334 r 14 - 1793232 r 13 + 1839416 r 12 - 2003712 r 11 + 2286224 r 10 - 2421504 r 9 + 3095088 r 8 - 
4428288 r + 5889152 r 6 - 6093312 r 5 + 4557056 r 4 - 2340864 r 3 + 774144 r 2 - 147456 r + 12288) 

(-12r + 7r 2 +4) 2 ]/ [8398080 r 10 ] 

where A(^>(A,Ni,Qi,L 3 ,N 3 ,N 2 )) = [\/3(-12xV - 24 x - 24 r - 12r 3 x + 9r 4 x 2 + 4y 2 - 8 \/3r y + 6 r 4 V^y x + 8 yfiy + 
+ 24 r x + 12 x 2 - 8 %/3y x + 4 r 2 y 2 - 4 r 3 \/3y + 3 r 4 y 2 + 4 V3r 2 y + 12 r 2 x + 12)] j [24 r 2 j . 



+ 11 ) 1(^)3 «W» - 



7 + 



> r- + : 



Case 13: 

x 3 } c NZ, E (Xi) u rKXj), Xi e 

+ / / . . + / / IFrnTw ^ = 



P({X 2 ,X 3 } C iV5js(Xi) U n(Xi), G T.) = 



where 

12r 2 + 24rx + 12x 2 - 
Case 14 



i ^(TO)) 3 

[l96608 - 3538944 r + 8927944704 r 7 - 1883996112 r 12 - 9492593152 r 8 - 146866176 r 3 + 29196288 r 2 + 220250112 r 4 - 
4486594560 r 6 + 213597 r 20 - 25 9 2 5 9 04 r 16 + 69124752 r 17 - 10683306 r 18 + 864387072 r 5 + 5220357120 r 9 - 

1081136256 r 10 + 602097408 r 11 + 2223664128 r 13 - 1509638512 r 14 + 716568768 r 15 ] j [l6796160 r 8 J 

N 1 ,Q 1 ,L 3 ,N3,N 2 )) = [V3(-12a; 2 r- 2 -24x-24r- 12r 3 a; + 9r 4 a; 2 +4i/ 2 -8\/3ry + 6r 4 % /3?/x + 8 % /3y + 
- 12 x 2 - 8 \/% x + 4 rV - 4 r 3 ^3y + 3 r 4 y 2 + 4 ^/lr 2 y + 12 r 2 x + 12)] j [24 r 2 ] . 

P({X 2 ,X 3 } C ^(XO U ri(JSTi), Xi G T.) = 

- [(r - 1) (-16384 + 278528 r + 215136 r 28 + 40 1 76 r 26 + 215136 r 29 - 3381264 r 27 - 2301952 r 2 - 99212040 r 23 - 

25050384 r 24 - 312101312 r 8 - 7215869272 r 13 - 147586784 r 9 - 42770432 r 4 + 12591104 r 3 + 114049024 r S + 345810944 r 7 + 

55914462 r 21 - 2082969096 r 17 + 43443459 r 18 + 826941555 r 19 - 641846754 r 20 + 20 9 9 3 06 16 r 22 - 23 2 9 6 3 72 r 6 + 
1311322268 r 10 - 3191747236 r 11 + 5434516904 r 12 + 7756861008 r 14 - 6865898928 r 15 + 4727296416 r 16 + 26115696 r 25 )] / 

[466560 (2r 2 + 1 n5 



\r y + 6 r 2 + 6 r x + 3 x 2 - 2 \/3y - 

■yj. 



where A(^(A, N 1 ,Q 1 ,L 3 , L A , Q 2 , N 2 )) = [\/3(-3 xV - 6 a; - 12 r - 6 r 3 x + 3 r 4 x 2 + 2 \/3i 

2 ^3r 2 y a; + 2 r 4 ^y a; + 2 \/3r 2 y - r V + 5 y 2 - 2 r 3 ^y + r A y 2 - 2 x + 6 + 6 r 2 x)] / 
Case 15: 

[(63855 r 10 - 498960 r 9 + 1650060 r 8 - 3036960 r + 3703292 r 6 - 3657696 r 5 + 3268368 r 4 - 2419200 r 3 + 1550448 r 2 - 

725760 r + 155520) (-6 + 5 r) 2 ] / [4199040 r 

where A(^(A, Ni, Qi, G 3 , M 2 , N 3 , N 2 )) = [ % /3(4r?/ 2 + 12a; + As/Zry + 9r - 4^32/ + 12ra: - 12 + 9 xV + 6 VE; 
12 rx 2 -4V3r 2 y- 12 r 2 x + 3 r 3 y 2 )] / [24 r] . 



\r 3 yx — 



Case 16: 



P({a 2 , A3} C JVp B (Xi) Uri(Ai), Ai £ T„) = / 17^7^ Vi3 dydx = 

Js 14 Jt 12 (x) Ay 1 \yvr 

- [(293 r 16 + 2344 r 15 +4662 r 14 - 9088 r 13 - 32320 r 12 + 42976 r 11 + 175408 r w - 119680 r 9 - 544144 r 8 + 372352 r 7 + 
1216512 r 6 - 882688 r 5 - 1564672 r 4 + 1373184 r 3 + 924672 r 2 - 1314816 r + 380928) 

(-2 + r) (r 2 + 2r-4) 2 ]/ [23040 (r + 2) 5 r 4 ] 

where A(^*(A, JVi, Qi,L 3 , N 3 , N 2 )) = [^(-12 xV - 24 x - 24 r - 12r 3 x + 9r 4 x 2 + 4y 2 - 8 V%r y + 6 r^^Jlyx + 8 V^y - 
12 r 2 + 24 r x + 12 x 2 - 8 \/3y x + 4 r V - 4 r 3 V3y + 3 r 4 y 2 + 4 \/3r 2 ?/ + 12 r 2 x + 12)] j [24 r 2 ] . 



Case 17: 

p({x 2 ,x 3 } c Jv; B (x,)ur;(ii), x 1 e t, 

r< " mW rl/2 ^ A(&>(A,N 1 ,Q 1 ,L 3 ,L 4 ,Q 2 ,N 2 )) 2 



/ / + / 

J sn J r i2 (x) J si 



A(T(y 3 )y 



[(6723 r 20 + 739 53 r 19 + 213678 r 18 - 433512 r 17 - 2873232 r 16 + 627264 r 15 + 20218896 r 14 + 5675184 r 13 - 97577924 r 12 - 
39916108 r 11 + 343932568 r 10 + 108508576 r 9 - 906967296 r 8 - 96480192 r + 1702951296 r 6 - 293251072 r 5 - 1994987520 r 4 + 

981590016 r 3 + 1118830592 r 2 - 1135919104 r + 287604736) (r - 1)] j [466560 r 4 (r + 2) 5 ] 

where A(0>(A, N x , Qi, L 3 , N 3 , N 2 )) = ^/3(-3x 2 r 2 - 6x - 12 r - 6r 3 a; + 3r 4 a; 2 + Isflry + 6r 2 + 6ri + 3i 2 - 2>/3j/ - 
2 ^3r 2 y x + 2 r*V3y x + 2 V3r 2 y - r 2 y 2 + 5 y 2 - 2 r 3 ^3y + r 4 y 2 - 2 v 7 ^ x + 6 + 6 r 2 x) J / \& r 2 j . 

Adding up the P({X2, A3} C Ap B (Ai) U r^(Ai), Ai £ T a ) values in the 17 possible cases above, and multiplying by 
6 we get, for r G [6/5, \/E — l) , 

^or(r) = - [-413208 r + 3070468 r 2 - 74801558 r 8 + 75243552 r 13 - 4883958 r 9 + 14541630 r 4 + 28880 - 11254002 r 3 - 
3667716 r 5 + 64360782 r + 13122 r 21 - 3300900 r 17 + 156014 r 18 - 175011 r 19 + 62825 r 20 + 1458 r 22 - 198 1 2000 r 6 + 
99831906 r w - 120628524 r 11 + 33155180 r 12 - 67685050 r 14 + 5055135 r 15 + 11053023 r 16 ] / [ll6640 r 6 (r 2 + l) 



(2r 2 + l) (r + 2) 3 (r + l) 3 ]. 



The v or (r) values for the other intervals can be calculated similarly. 



Appendix 4: Proof of Corollary I5.lt 

Recall that S^ nd (r) = Pj nd (r) is the relative edge density of the AND- underlying graph for the multiple triangle 
case. Then the expectation of S^ nd (r) is 

e [srV)] - n J_ x) E E E K d ( r )] = E WW] = p(x 2 g AW*i) n rj(Xi)) = p and (r). 
n ^ > i< 3 

But, by definition of N r PE {-) and T\{-), if A\ and A~ 2 are in different triangles, then P(X 2 G A^A^nr^Ai)) = 
0. So by the law of total probability 

Pand(r) := P(X 2 eN^ E (X 1 )nT r 1 (X 1 )) 

= E p (* 2 e iVp E (^i)nr[(x 1 )|{x 1 ,x 2 } ctap^Xx.x,} ciy 

i=l 

■/„« 

= E^ d M ^({^i, X 2 } C Ti) (since P(X 2 € A^ £ (*i) n rjpfi) | {X l ,X 2 } C T<) = p and (r)) 
i=i 

E ( E ffU ) (sincc p({Xl '* 2} c r <> = ( E ^ (r ^ ) ) = a-M (l>i ' 

where Pand(^) is given by Equation (|T^|) . 

Likewise, we get p r(f) = Por(?*) w f ) where p OT {r) is given by Equation (|T3l) . 

Furthermore, the asymptotic variance is 

£ and (r) = E [Kf(r)htf(r)] - E [/^ d (r)] E [/^ d (r)] = P({X 2 , X 3 } C JV^pfO n rftXi)) - (p a „d(r)) 2 . 
Then for J rn > 1, we have 



2 

i=l 



p({x 2l x 3 } c ^(xonru^)) = E p (^,x 3 } c Ar^(x 1 )nri(A' 1 )|{x 1 ,x 2j X3} c T i )P({x 1 ,x 2 ,x 3 } c r<) 

i=l 

= P({x 2 , x 3 } c ^(Xr) n rj(*i) I {Xi, x 2 , x 3 } c r e ) ( E ^ 3 J ■ 



Hence, 



v,i=i 



?a„d(r) = p({x 2 ,x 3 } c N£ E (X!) n r;(JTi) | {*i,x 2 ,* 3 } c r e ) (XX ) ~ (p a „d(r)f 



Likewise, we get i/ or (r) = i> or (r) (e/=i w ?) + Por(r) 2 f Ef=i ^ 3 ~ (E/=i w ? 



So conditional on y m , if i? an dM > then -y/n (S^ nd (?') — p a nd(r)) — ► N (0, ? a nd(r)). A similar result holds 
for the OR-underlying version. ■ 

Appendix 5: Proof of Theorem 15.21 

Recall that pj " d „ (r) is the version II of the relative edge density of the AND- underlying graph for the multiple 
triangle case. Then the expectation of pf} d n {r) is 



2n t 



since by ((3]) we have 



riiirii - 1) 



and / 
kl \ 



P T(r) =Pa„d(r) 



E[^ d (r)] =p a „d(r 



k<l 



where p a nd(r) is given by Equation (|12[) . Likewise, we get p OI (r) = p OI (r) where p or (r) is given by Equation (fT3]) . 
Next, 



Var [p?? d (r)] = £ 



^ 2 K - i) s 

4ii? 



Var 



since /?f" d (?") and pff d (V) arc independent for k ^ I. Then by Q we have 



Var 



Tli (Hi - 1) 



4 fn ■ — 2) 

Var [^ d (r)] + \ 1 _ Cov [^ d (r), /i^ d (r)] 



So, 



Var [^ d „(r)] = £ Var [/^ d (r)] + £ 



m {ni - l) (m - 2) 



2^ 



Cov[/^ d (r),^£ d (r)]. 



i=i 



Here "' Var [ /l i2 d ( r )] = ^7 Var [^M] ■ Thcn for lar S c m and 

-Var [/^ d (r)] « — ^ -Var [^ d (r)] 

«t n 2 Ei=i < 



since ^\ = J~)j—i "' ^ and rii/n Wi as rij , n — > 00 . Similarly, for large n, and n, 



E 



rii (m - 1) (rij - 2) 



Cov [AS d (r),feS d ( 



W E' 



iJ =1 



Cov[^ d (r),^ d (r)] 



Hence, conditional on y m , sjn (/of/ d n (r) -pand(r)) — > A" (0, 4 i/ and (Y)) provided that £ an d(V) > where 
Pand(V) = Pand(r) and u an d(r) = iw(r) (j2i=i w i) / (j2i=i w f) ■ A similar result holds for the OR-undcrlying 



version. 




Figure 15: An illustration ol the vertices for possible types of N PE (xi) D T^(xi) for x\ € T s . 




Figure 16: Prototype regions Ri for various types of Np E (x\) n rj(xi) and the corresponding points whose 
^-coordinates are Sk values. 



